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FOREWORD 


This research was conducted by The Pennsylvania State. University 
for the National Aeronautics and Space Administration, George C. Marshall 
Space Flight Center, Huntsville, Alabama . Dr. George H. Fichtl of the 
Aerospace Environment Division, Aero-Astrodynamics Laboratory, was the 
Scientific Monitor. Professor John A. Dutton, Department of Meteorology, 
The Pennsylvania State University, was the Principal Investigator. The 
support for this research was provided by Mr. John Enders of the Aero- 
nautical Operating Systems Division, Office of Advanced Research and 
Technology, NASA Headquarters. 

The research reported in this document was motivated by the need 
to indicate critical atmospheric parameters in relation to associated 
critical aircraft response characteristics for the safe design and opera- 
tion of aeronautical systems. This res^earch is part of an investigation 
of the effect of severe atmospheric forcing on aeronautical systems. 

The report is a Master of Science Thesis in Aerospace Engineering prepared 
by Mr. Hillard in collaboration with and under the direction of 
Professor Maurice M. Sevik, Professor of Aerospace Engineering, The 
Pennsylvania State University. 

A previous report presented results of an empirical study of 
thunderstorm gust fronts, and a subsequent report will present results 
from a numerical model of gust front structure and evolution. 
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1.0 ISTRODUCTION 


1.1 Previous Investigations 

The local fluctuations of air velocity encountered by a craft 
moving through atmospheric tvirbulence are random fiinctions of position 
and time which are computable only in a statistical sense. They, 
therefore, produce aerodynamic forces and moments which are likewise 
expressible only in a statistical sense. If the velocity field can 
be approximated by one for which the turbulence is statistically 
stationary and ergodic, mathematical techniques have been developed 
which may be used to find the rigid body and elastic response of the 
aircraft [l]. 

Several investigators have treated the problem of gust en- 
counter. The first significant contributions were made by Wagner [2] 
1925 , who considered the problem of a step change in angle of attack. 
Kussner [3] 1932, considered lift generation from a sharp edge gust. 

A fundamental method was presented by von Karman [U] 1934 for the 
treatment of arbitrary airfoil motion. In his method, von Karman 
determined the system of impiilses which are applied by the wing to 
the fluid. He accomplished this by investigating the three dimensional 
vortex wake phenomenon. Sears [ 5 ] 19^0, introduced a Laplace trans- 
formation on the time variable which greatly simplified the numerical 
integration which was necesssury in the use of von Karman 's method. 

Sears also calculated a function which predicts the lift due to a 
sinusoidal gust. The general case of an inclined sinusoidal gust has 
been treated by Filotas [6] 196? who obtained a closed fonn expression 



for the loading on a two dimensional airfoil. Horlock [j] 19^8 treated 
the component of gust which acts parallel to the wing chord and 
developed a function which yields the contribution to lift of this 
component. Naumann and Yeh [8] 1972 have added the effect of a cambered 
airfoil to Horlock 's analysis. Sear's function has been used in con- 
junction with two dimensional lifting surface theory by Giesing [9] 

1969 for the treatment of harmonic gusts. 

The finite wing has been treated by several investigators. 

Jones [10] 19^0 calculated the derivative coefficients for airfoils of 
finite span and also for airfoil-aileron-tab combinations. Reissner [ll] 
19^3 studied the effect of finite span upon the airloads on oscillating 
wings for subsonic compressible flow. In an unpublished paper by 
Graham [12], he succeeds in obtaining a numerical solution for the 
load on a thin rectangular wing in an arbitrary incompressible flow 
field. 

Random atmospheric forcing is generally treated using power 
spectral techniques. Diederich and Eggleston [13] 1956 calculated the 
power spectra of moments on a wing in random turbulence neglecting 
the effect of side gusts and spanwise variations in gust intensity. 
Diederich and Drischler [l4»l6] 1957 included spanwise gust variations in 
their investigation of sharp edged traveling gusts and the lift occurring 
from penetration of these gusts. Houbolt [15] 1970 has added important 
procedural methods for establishing gust load levels appropriate for use 
in aircraft structural design. A method of establishing load exceedance 
curves for structural fatigue analysis was presented in this study. 
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The ansJ-yses due to Diederich and Efegleston [12], and Diederich 
and Drischler [13], contain several fundamental assumptions. It is 
assumed that the encoxintered t\irbulence is homogeneous and isotropic. 
It is fiurther assumed that the turhiolence displays gust intensity 
invariance along the flight path until the craft has traversed a 
given region of tiirhulence. Two dimensional, vinsteady aerodynamics 
for incompressible flow was used in these studies. Although the above 
assumptions are compatible to certain cases, their restrictions limit 
the generality of the results obtained. 

Houbolt [ 15 ] 1970 , used power spectral techniques to treat the 
problem of vertical gusts that are uniform in the spanwise direction. 
The atmospheric turbulence environment is modelled in terms of a 
spectrum shape. The parameters of this model are the r.m.s. gust 
intensity, turbiilence scale length, and a parameter which allows for 
the flight time spent in the tiu-bulence. The response is given in the 
frequency domain. Gust load levels are established for aircraft 
structural design. Houbolt 's work provides systanatic procedxires for 
structural response to random conditions with no more difficulty than 
is usually associated with a discrete gust analysis. 

1.2 Investigation of Severe Atmospheric Forcing 

Severe atmopsheric phenomena affect several aspects of aircraft 
design. Adequate protection to systems must include the effects of 
severe conditions on returning space vehicles, vertical take-off 
aircraft, automatic control systems, and their respective structiiral 
integrities . 
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Statistical techniques now being utilized to investigate the 
response of a vehicle to atmospheric forcing make two basic assumptions; 
(a) the system is linear and (b) the input is realized from a stationary 
Gaussian process. These assumptions are not always true. 

The investigation of atmospheric forcing necessarily includes 
the basic properties of the craft's response. The correlation of the 
response to the forcing function provide valuable insight in two 
important areas . First such analysis provides information to the systems 
designer, and secondly it provides the meteorologist with clues as to 
what turbulence parameters most strongly affect such systems. It is 
with these two factors in mind that this study is undertaken. 

1.3 Statement of the Problem 

Designing for unusual atmospheric conditions requires the 
determination of structural stress levels, rigid body motions, and 
elastic responses of the craft. The aerodynamic forces emanating 
from random forcing are comprised of a component due to the external 
gust input and a component due to the craft's rigid and elastic 
responses to such input. If these contributions are calculable, then 
the responses can be found. 

The specific problem to be considered in this study is that of an 
aircraft in steady level flight encountering a spanwise gust wave. 

Figure 1.1. The previous investigations of Diederich and Eggleston [13], 
Diederich and Drischler [lU], and Houbolt ]15] have dealt with the 
correlation functions and power spectra of the forces and moments which 
occur to a craft in a turbulent environment. This study treats those 
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aircraft parameters and meteorological p8u:ameters which are most critical 
to aircraft response. Results are presented which illustrate the depen- 
dence of these critical parameters with respect to Mach number , relative 
mass, and nat\iral frequency. In order to reduce the mathematical com- 
plexity of the problem without compromising the significance of the 
results and conclusions, the craft's wing shall be unswept and untapered. 



2.0 THE RIGID BODY EQUATIONS OF MOTION 


2.1 Review of the Eqimtions of Motion 

The force and moment acting on and about the center of mass of an 

aircraft are directly derivable from Newton's second law. The center 

of mass has the linear velocity components U, V, and W, and angular 

momentiim components h , h , and h . Their time derivatives are given by 

X y z 

dVc dVc 

dt “ + 0) X Vc 

dh 6h 

at ■ at * Si " ii 

where the coordinate system is attached to the craft and in motion with 

it. From expressions (2.1) and (2.2), the equations describing the 

* 

plunging ancT pitching motions of an aircraft become: 

(2.3) F = m (W + PV - QU) 
z 

(2.1») M = h + Rh - Ph 
y X z 

where P, Q, and R are the angular velocity components. 

The forces acting on an aircraft are gravitational and aero- 
dynamic in nature. The usueG. assumption may be made that the 
disturbances are small deviations from the reference flight state. 

This assiimption gives a reasonable approximation for disturbances to 
9 less than •rr/6. Generally flight of considerable turbulence occurs 
with quite sTn«-H linear and angular magnitudes. The variables may 
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*be assumed to be the sum of their reference flight value and a small 
distixrbance value, thus enabling the linearization of the equations 
of motion: 

(2.5) AZ = m ^ - mu ^ 

dt^ ° 

(2.6) AM = B ^ 


The equations may be expressed in nondimensional form by dividing 

o p _ 

( 2 . 5 ) and ( 2 . 6 ) respectively by 1/2 pu S and 1/2 pu Sc. The resulting 

o o 

expressions are given by, 


(2.T) AZ = yD (Dh -6) 

(2.8) AM = 2i_D^e 

B 


2.2 Relations of Unsteady Lift and Moment 

The unsteady lift and moment forcing functions must be expressed 
with respect to a set of coordinate axes fixed in the aircraft as 
follows : 

( 2 . 9 ) = - TTpbV I - - (|- + a) a} 

/N 

(2.10) » TTpbV (|- + S) ^ 

+ (i + a)^ L ] a} 

^ n 

These relations then give the component of unsteady force and moment 
due to the aircraft’s response. The lift generated on an airfoil due 
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to a sinusoidally varying gust may be expressed in terms of Sears 
function $(k)i it is located at the quarter chord point. The lift 
per unit span is given by 

(2.11) = 2Trpu^ba_$(k)e^‘*’^ 

o u 


In equations (2.9) and (2.10), the lift and moment are applied at 
and about the elastic axis. The aerodynamic coefficients and Sear's 
function appearing in these equations have the numerical, values shown 
in Smilg and Wassermann. 

In order to obtain the total force and moment acting on the craft , 
the effects of the above stated forces and moments must be summed over 
the wing and tail. The total force AZ and moment AM may be expressed as 


(2.12) AZ = 2 


r-M ,.G _M 

[L + L + L . + L. Jdy 
w w t t 


(2.13) AM = 2 


[M„ . . X^L« - x;l° . 


The subscripts on the moment and lift terms indicate the wing 
or tail contribution. The superscript indicates the source of the 

/S 

term to be either gust or motion. The terms indicate geometrical 
distances from the center of gravity to the points of application of the 
forces. These distances are shown in Figure 2.1. 

Substituting equations (2.12) and (2.13) in nondimens ional form 
into equations (2.7) and (2.8), the rigid body equations of motion are 
obtained. In Chapter U.O of this study it will be shown that the tail 
frequency parameter k. is equal to 0.6k for a large class of subsonic 
aircraft. At this time the relation will be used without proof. 
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Figure 2.1 - Location of Geometrical Distances 


Note should also be made that at any given instant of time, the 
tail forces will have a phase difference with respect to the wing 
forces. A simplification to this situation will also be given in 
Section 2.6. 

2.3 Solutions for the Unknown Normal Coordinates 

The general form of the solution of the equations of motion may 
be expressed as the sum of the products of the mode shapes with cor- 
responding generalized coordinates as follows. 

(2.1U) h = Z:(|)j(y)Cj(t) 

Since varies sinusoidadly with time, it may be expressed 
J 

in terms of the mechanical admittance of the system as follows : 

(2.15) 5j(t) = Hj(i0))e^“^ 

Therefore, eqviation (2.lU) becomes 

( 2 . 16 ) h = 2(j)j(y)Hj(ia))e^“* 

J 

The mechanical admittance functions occurring in the rigid 
body equations of motion contains four terms, namely: plunging H^(io)), 

elastic bending H 2 (iw), pitching T^(io)), smd torsion T 2 (iw). 

The coordinates pertinent to these equations and their 
derivatives may then be expressed as 
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(2.17) = Z(J.j(y)Hj(i(u)e 


io)t 


(2.18) h, = 


2<(>j (y)Hj (iw)6^je^‘^^ + X2li|»j(y)Tj(y)Tj(ia))6j^je 

J J 


iwt 


(2.19) a = 


w 


+ liwi 26 (y)H,(ico)e 
o 0 j 


Icot 


+ X Eij; (y)T (iu)e^“^ 

^ Uq j J J 


( 2 . 20 ) = 


(—) + 


j. ( ) T"j. / _\tt ^ _• . . \ r _iWt 


“o “o 


^ J 


(2.21) 0 = ZiPj(y)Tj(io))6^je 


iwt 


(2.22) 0 = (ia))^2i|;j(y)Tj(iu))6^je^‘^^ 

J 


The term 0 represents the rigid body angular displacement coordinate. 

The term a represents the total angle of attack. This includes in 

addition to all response coordinates, the effect of the vertical gust 

and plunging velocity on the angle of attack. The term equals 1 

when j=l and zero otherwise and signifies that only rigid modes are 

included. It should also be noted that w * includes the phase 

g 

difference which the gust velocity at the tail has with respect to 
that acting on the wing. While the equations themselves describe 
the rigid body behavior of the craft, the forcing functions contain 
elasticity and torsional contributions to this behavior. 
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2.k Rigid Body Ecruatioits of tfotion 


The equations of motion for the rigid body response of the aircraft 
may now be stated. The following substitutions for terms containing 
aerodynamic coefficients are made: 

( 2 . 23 ) L^(k) = [L^ - (|-+ a) 

(2.21+) LgOt) = - (|-+ 5) 

(2.25) L^(k) = [M^ - (|•.+ i) (L^ + 

( 2 . 26 ) L^(k^) = [L^ - (| . 

(2.27) Lj(k^) “ ^ 

(2.28) Lg(k^) = [M^ - (i . S^) (L^ . M^) + (| * S^)^ 

The equation for rigid body plunging is given by 


o 


<1)t (y) 

H^(ik) {l+(ik)^ 

[y - irdij^ + (ik) L^) 

il 

p 

A 4 > (y) 

O 

- TTb^ (L, 

* (iS) i 

dy]} 


o 

fH- <|)p(y) 

+ H2(ik) { 

-UiT (ik)^ (Lj^ + (ik) : 

j i 

^ O 


p 

iI't (y) 

+ T^(ik) 

{-2(ik) [y + 2 tt (ik) 

X3L1 , dy 


5 - (y) 

+ 2irb‘^(^) (I^ + b(^) Lj^) J <3y]} 

+ T 2 (ik) {UiT(i^)^ XgL^ I — dy} 
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= 2tt (^) {[(ik)^ + 2$(k)] 

o 


I 


+ [b^(i^)^ Lj^ + 2$(k^)] ^ 


Similarly the equation for rigid body pitching motion may be written 
, as follows : 


(2.30) H^(ik) {2Tr(i^)^ [L^ + (ik) + 2 X 2 \ 


+ 2(ik) X2L^] 


^ 4>3^(y) 
£ 


dy 


+ 2iT(ik)^ b^ [b + b^ (ik) Lg + 2X^ 


+ 2(i^) b 


\ ^i(y) 


dy} 


+ H 2 (ik) {2ir(ik)^ [L^ + (ik) + 2 X 2 L^ + 2(ilt) X 2 L^] 


<|)o(y) 


dy} 


+ T^(ik) {2(i^)^ [Uig + Tr(ik) {(X 2 L 2 + ^^2^1^ 




dy 


(b^ X^L^ + 2b^ X^Lj^) 




dy]} 


+ T2(ik) {2iT(d^)^ X2 (L^ + 2X2L^) 




dy} = 


- TT(^)e^'^ [4X^b$(k^) + b^(ik)^ Lg + 2b^(ik)^ X^L^] 


t ^ 
£ 


w 

-TT(ik)2 (^) [L^ + 2X2L^] 

o 


£ 
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The equations as derived provide two relations necessary for 
computation of the longitudinal rigid response of a craft. The inclusion 
of bending and torsion modes necessitate the introduction of additional 
equations, the number of which will correspond to the number of elastic 
modes considered. 

2.5 Effect of Downwash 

The aerodynamic characteristics of the component parts of an air- 
craft are affected by each other. These interferences occur because 
of interactions between such components as the wing and fuselage, tail 
and fuselage, or the wing and the tail. 

The most important effect is the downwash contribution that is 
due to the unsteady forces acting bn the wing of the aircraft. Since 
the unsteady lift and moment expressions, (2.8) and (2.9), are 
dependent upon the instantaneous euigle of attack of the tail, it is 
necessary to decresise this angle by the amount of downwash angle e. 

A general expression may be written for the unsteady moment acting 
on the craft, namely 

(2.31) = (Xj^L^ - XgL^) cos - (X^L° + X^L^) cos (a^ - e) 


+ M + 
w t 


The angles and (a^ - e) may be assumed small. Equation (2.31) 


may therefore be written in nondimens ional form as 


(2.32) C = (Xj^C® - X - X C° - X 3 (^ + = C ) 

u t t t 
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The terms of equation (2.32) which must he investigated are 
and C„ , and the significance of the downwash must he determined. 

The time dependent downwash angle e may he approximated hy the 
expression, 


(2.33) 


We substitute in equation (2.32) the expressions for lift and 
moment, namely equations (2.9), (2.10), and (2.11) and obtain the 
following expression for the unsteady moment coefficient: 


(2.3H) 


= H^(ik) ^ h * ^2 


1 f^<t>i(y) o ^ 

- |(ik) L^] J ^4“ f^^^3 ^h ^ ^5 


+ (ik)h (1 - 1^) (Lj^ + - Lg)] — j— 

^ o 


p - 1 '*’2^^^ 

+ H2(ik) {iTk [(2X^1^ + Lg) - (ik)(2L^ 2 J ~~l — 


+ T^(ik) {-TTk^(ik) Xg (2L 




a 4^i(y) 


+ b^X^k^ f ^ i ^"6^ 

"t 




r^- i|^o(y) 


+ Tg(ik) {-TTk^(ik)Xg (2L^ I J ~ 


+ {Tr[-k^ (^) L, - ^kV (-2-) + 2X,W $(k)] f- 


'u ' 1 2 3 u 

o o 
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[TTk^ 


-3 


^4 

o 


w i-4 , 2 


(r:^) LJ 




(1 


9£ v f^t ^ 

3a^ 1 I 

■'o 


As 

- 2ttX b 

X i5 





This equation shows that the downwash effects appear only in the 
rigid body modes. A conservative estimate of the effect of downwash 
may be established by considering the relative change to equation 
(2.34) upon deletion of 3e/3a. The contribution due to downwash effects 
will be of the order of less than ten percent for the unsteady moment. 
For conventional aircraft, a typiceil value of de/da is 0.40. 

Similsurly for the downwash contribution due to lift, the lift on 
the aircraft may be expressed by: 


(2.35) + lJ + l” + l' 


w 


which in nondimen&ional form hecomes equal to 


(2 


.36) = H^(ik) {TTk^ [(21^ - (ik)L^^) j 


^ 4>3_(y) 


dy 


+ (2b^ + b^(ik) Lj^ (l - 3e/3a)) 




dy]> 


o 

p < t > p ( y ) 

+ HgUk) {TTk'" [21^ - (ik) L^] dy} 

O 


/s. fS, \j)^ (y ) __ ^ 

+ Tj^(ik) {-TTk^ [(ik) ^ — dy + (b"^(j^) L^X^ 

J r \ 


+ 2b^ 


^3\> - to ’ 
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+ T^Cik) {-trk^(ik)X2L^ 


A 'i'o(y) 


dy} = 


w 

{TTbe^*^ [2w^H\) + k^% (^) (1 “If)] 

o 


t ^ 
Z 


"W 

+ tt[2w $(k) - k^ (— ) L ] 
u U 1 

O 


I ^ 


Evaluation of equation (2.36) for the effect of downwa'sh yields 
a contriuution of not more than twelve percent of the unsteady lift 
coefficient. It may, therefore, he concluded that the effects of 
downwash may he safely neglected in this analysis. 

2.6 Wing to Tail Phase Difference 

Due to their spatial separation a phase difference $ will occur 

between the wing and tail. For a given atmospheric disturbance of 

wavelength 3, the time required to traverse a cycle is B/u^. Similarly, 

the distance from wing to tail is approximately and the time 

required to traverse this distance X^/u . In terms of the frequency, 

o o 

the phase difference may he expressed as 

(JL)X 

(2.37) ^=— 

o 

which may he written in the form, 

( 2 . 38 ) $ = 2)s.X^ 

As is shown in Chapter U.O, the nondimens ional distance is 
a constant over the entire range of aircraft considered. The reduced 
frequency k therefore dictates the value of the phase difference. 
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3.0 EIASTIC EQUATIOMS OF MOTION 


3 . 1 Introduct ion 

An aircraft is a single elastic body. The interactions of the 
elastic components present many complications in the analysis of such 
a structure. It is therefore analytically convenient to consider a 
craft as either rigid or semi rigid. In many cases this procedure 
will yield a reasonable engineering approximation. 

For this analysis the craft will be asstmed to be rigid except 
for the inclusion of wing elastic bending and torsion in their 
fundamental modes. 


3.2 Lagrange^ s Formulation of the Equations of Motion for an 
Elastic Wing 

Assuming the actual motion of a wing to be comprised of a 
combination of fundamental wing bending and wing torsion, the system 
may be expressed as an equivalent system consisting of an airfoil 
section of unit span restrained by springs against bending and 
torsional motion [1T]- 

Using the Lagrangian approach, the equations of motion are 
given by 

(3.1) ME + S^oi + Mw^h = 


(3.2) 


S h + I a + I w a = 
a a a a a 


where the generalized force and moment and represent the unsteady 
forcing to the wing per unit span. The generalized forcing terms may 
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be expressed using the quantities which have been previously defined in 
Chapter 2.0 and are restated here for convenience: 


(2.9) = - TTpbV ■ ^2 ^h^ 

/N 

(2.10) = TrpbV I - - ^2 “ ^n^ 

+ (| + Sj L ] a} 
d — n 

(2.11) = 2TTpu^ba„v(K)e^'^^ 

o u 

The generalized force and moment are, therefore, of the form 
(3.3) S. ' 

(3.M = L° . My 

These equations may be nondimensionalized in a manner analogous 
to the rigid body equations of Chapter 3.0. One unique distinction 
should be made, however. The left hand side of equations (3.3) 
and (3.U) describe only the elastic bending and torsional modes and 
does not contain any rigid body considerations . The rigid body mode 
has previously been designated as the first mode thus necessitating 
the introduction of a Kroenecker delta which was written as • For 
the elastic modes, the subscript i ranges from 2 to n, depending upon the 
nmber of modes considered. 

3.3 Solution for the Equations 

The elastic deformations are given by the sum of the products 
of the mode shapes and the corresponding normal coordinates as follows : 
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(3.5) 


\ ~ Hj(iw)e 


Ktft 


(3.6) 


(3.7) 


J^<l>.(y) Hj(iu))e^‘^ 


" "o % P 


+ X. 


^2 "T^ T (iw)e 

o J J 


i(i)t 


hg = Z<{'j(y) H^(iw)Sj^ e 

<3 J 


ili)t 


(3.8) 


h = (^)2 z<j, (y) H,(iw)5. 
J 


(3.9) a = Ei|» (y) T (iw)5 e 

J j ij 


io)t 


(3.10) a = (io))^ 24>.(y) T (iw)5. 

P •’ 


where h is the wing pliuaging coordinate and h^ is the translational 

coordinate due to wing bending.. The term is the wing angle of attack, 

and ttg the contribution to angle of attack due to wing torsional motion. 

The bending and torsionsQ. mode shapes are represented as <t*.(y) and ij^.(y). 

J J 

Hj(io)) and Tj(iu)) are the admittance functions for the translational 
modes and torsional modes respectively. 

These equations must next be integrated over the semi span length. 

It is advantageous - prior to the integration - to multiply the equations 
l^y <l>j^(y) for the bending equation and ilf^(y) for the torsional equation. 
The orthogonality of the modes then yields : 
fA 


(3.11) [ m ((>.(y) <(>.(y) dy = M6. 

Jo ^ 
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(3.12) 


i 4>Ay) (y) dy = i6. 

>o ^ 

where 6^j has its usual definition. 

3.4 Equation for Elastic Bending 

Substituting equations (3.3) through (3.10) into (3-1) and (3.2), 
we obtain the equation for elastic bending: 

(3.13) H^(ik) {-2TT(ik)^ a^^ [l^ + (ik)L^]} 

+ H^(ik) {2a [U(ik)^ m - 4(ik) )^^ in 

^ ij — — h 

-7r(ik)^ - Tr(ik)^L^]} 

+ T^(ik) {-2 (ik)^ ^ 

+ T^(ik) {2(ilt)^ [4S - ir(dlc) 

w 

= TT {(^) g. [(ik)^L^ + 2$(k)]} 
o 

3*5 Equation for Wing Torsion 

Similarly the equation of torsion may be expressed as 

(3.14) H^(ik) {7T(ik)^S^j [L^ + (ik)!^]} 

+ H2(ik) {S^j(ik)^ [8S^ + ttL^ + Ti(ik)L2]} 

+ (ik) {TT(ik)3 X L- t. .} 

1 — 2 3 ij 
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The form of equations (3.13) and (3.lU) indicates that they comprise 
a set of equations, the number of which corresponds to the nimiber of 
modes which are to be included. The set of equations established now 
enables the solution of the functions Hj(ik) and Tj(ik), for any desired 
number of elastic and torsional modes. 



h,0 AIRCRAFT PARAMETER CHARACTERISTICS 


h . 1 Introduction 

Examination of the equations of motion indicate that there exists 
a large number of variables. It would be computationally advantageous 
to reduce the number of parameters to a minimum. One would then be able 
to study how the various dimensionless parameters relate to each other 
for a large set of subsonic aircraft. 

The equations of motion derived for an aircraft in Chapters 2.0 and 
3.0 show dependence on aerodynamic coefficients and also on parameters 
which are directly attributable to the mass and geometrical characteristics 
of the craft. General observation of a particular class of aircraft 
indicate that they display similar basic characteristics. These similar- 
ities emanate from the function of the craft and the aerodynamic, control, 
and structural qualities which are necessary to accomplish this function. 
Since the parameters of a craft are relatively closely dependent upon 
its function, with the primary difference being the overall size of various 
craft within the^ class, one suspects that the various parameters might 
display some trends when considered over the entire range of that type 
of aircraft. Obviously the parameters should be plotted against a value 
indicating size. Gross take-off weight has been closen as most convenient. 

The utility of displaying such characteristics derive’s from the 
ability of a designer to interpolate or extrapolate the trend to obtain 
a close approximation of the value for any craft desired. This can then 
be used to obtain useful information for the design of new aircraft or as 
is used in this study, the computation of the response of the craft to 
some input condition. 
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It. 2 Classes of Parameters 


As an example, the class of subsonic transport jets shall be con- 
sidered. In general, the parameters displayed in the equations of motion 
fall into the following three distinct categories: 

(a) Parameters which are dependent upon mass or mass distribution 
of the craft. 

(b) Parameters which are dependent upon a geometrical length of the 
aircraft, but are not a function of spanwise location. 

(c) Parameters which are a nondimens ional geometrical length of 
the aircraft, but are a function of spanwise location. 

In the analysis as presented in this study, the airfoils were 
assumed to be unswept and untapered thus eliminating the third type of 
parameter. It should be noted, however, that an extensive analysis of such 
craft will contain this type of parameter. A discussion of such parameters 
appears in Section 5 ot this chapter. 

The class of parameters dependent on mass and mass distribution 
contain y the relative mass parameter, i„ the mass moment of inertia 
parameter, S which is the wing static moment term, m the nondimensional 
wing mass distribution parameter, and I the wing mass moment of inertia 
parameter. The nondimensional wing mass distribution terms are derived 
from the inclusion of elastic bending and torsional effects in the analysis. 
It is found that these quantities, due to their overall mass dependence, 
may be expressed as linear function of the relative mass parameter y. 

It is, therefore, necessary only to specify y in order that these 
quantities be defined. 
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The second class of parameters contain lengths or ratios of 


lengths in nondimensional form. These quantities are b, the ratio of 

/s XV xv 

tail to wing mean aerodynamic chord, and X^* X^* and X^^ which are 
the nondimensional distances from the aircraft’s center of gravity to the 
respective locations on the wing and tail where the unsteady aerodynamic 
forcing occurs. 

As previously indicated, parameters of the third class have not 

been used in this analysis with the exception of the mode shapes 

(j) (y) and 4^.(y) which occur because of elasticity considerations. 

J J 

4.3 Parameters of the First Type 

As previously mentioned, parameters of the first type may be 
defined by the specification of the relative mass parameter y, and 
the mass moment of inertia parameter about the spanwise axis i^, both 
representing nondimensional quantities. 

The relative mass parameter is defined by the relation 


W 


(U.l) y = ^^2 pg g - 


The mass moment of inertia parameter is similarly defined as 


ik.2) 




B 

3 

p Sc-^ 


When the parameters defined by equations (U.l) and {h.2) are 
plotted against gross weight, the results are shown in Figures U.l through 
4.5. Table 4.1 lists the crafts and their respective data which were 
used for parameter determination. Also shown in Figure 4.1 are the 
values of those wing parameters which can be seen to be linear functions 
of the relative mass parameter. 
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Figure U . 1+ - Wing Stat 
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Table k.l Pertinent Aircraft Data 
Aircraft Designation T.O. Weight 


1. 

Lear Jet 2k 

13,000 

2. 

Aero Jet Commander 

i6,8oo 

3. 

Dassault Fan Jet Falcon 

26,kp5 

U. 

Lockheed Jet Star 

1+0,921 

5. 

BAG One Eleven-200 

78,500 

6. 

Boeing 737-200 

93,500 

7. 

TU 13kA 

97 ,000 

8. 

Sud Aviation-210 Car. 

101,1+13 

9. 

Hawker-Siddeley Trident 

115,000 

10. 

Boeing 727-200 

160,000 

11. 

Boeing 720 

229,000 

12. 

Boeing 707-120B 

257,000 

13. 

Douglas DC-8-10 

273,000 

Ik. 

BAG VClO-1106 

312,000 

15. 

Douglas DC-8-50 

315,000 

l6. 

Lockheed Starlifter 

316,600 

17. 

Boeing 707-320B 

328,000 

18. 

Boeing 7^7 

680,000 

19. 

Lockheed C-5A 

728,000 


32 



Table U.l (Contixmed) 



2 

Wing Area ft 

Tail Area ft^ 

Wing M.A.' 

1. 

231.77 

n.a. 

7.00 

2. 

303.30 

n.a. 

7.67 

3 . 

U52.O 

n.a. 

9.33 

k . 

51+2.5 

n«a. 

10.10 

5 . 

980.0 

257.0 

12.00 

6 . 

980.0 

n.a. 

11.28 

7 . 

1,238.0 

n.a. 

19.90 

8 . 

1,579.0 

n.a. 

15.38 

9 . 

1 , 358.0 

310.0 

17.20 

10 . 

1,700.0 

376.0 

18.18 

11 . 

2 , 1 + 33.0 

500.0 

20.80 

12 . 

2,521.0 

51^5.0 

20.92 

13 . 

2 , 773.0 

n.a. 

22.20 

ll+. 

2,932.0 

638 . 0 

n.a. 

15 . 

2 , 773.0 

n.a. 

22.20 

16. 

3,228.0 

n.a. 

22.20 

17 . 

3 , 010.0 

625.0 

23.80 

18 . 

5 , 500.0 

1,1^70.0 

38.30 

19 . 

6 , 200.0 

n.a. 

33.20 
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Table U.l (Continued) 


Tail M.A.C. ft. Wing Semi -Span Tail Semi-Span 


1. 

3.50 

17.79 

7.33 

2. 

3.90 

21.65 

9.66 

3. 

3.63 

26.25 

10.66 

1+. 

6.72 

27.20 

12.38 

5. 

9.35 

44.25 

i4.75 

6. 

9. ho 

46.50 

18.00 

T. 

8.90 

47.57 

15.08 

8. 

9.20 

56.25 

17.38 

9. 

9.38 

44.92 

17.12 

10. 

11.35 

54.00 

17.88 

11. 

14.00 

65.42 

19.83 

12. 

13.50 

65.42 

21.67 

13. 

10.90 

71.20 

23.75 

lU. 

14.10 

73.08 

21.92 

15. 

10.90 

71.20 

23.75 

l6. 

n.a. 

79.96 

25.19 

IT. 

14.20 

72.88 

22.88 

l8. 

23.50 

97.83 

36.37 

19. 

14.15 

111.35 

34.35 


Table U.l (Continued) 


Tail to Wing 

Chord Ratio Wing Taper Ratio Rel. mass y 


1. 

0.500 

0.508 

566 

2. 

0.510 

0.333 

510 

3. 

0.602 

n.a. 

443 

U. 

0.665 

0.373 

529 

5. 

0.780 

0.323 

470 

6. 

0.830 

0.339 

590 

7. 

0.1^46 

0.221 

275 

8. 

0.598 

0.354 

295 

9. 

0.546 

0.268 

347 

10. 

0.625 

0.304 

365 

11. 

0.672 

0.333 

325 

12. 

0.645 

0.332 

343 

13. 

0.491 

0.230 

314 

lU. 

n.a. 

n.a. 

n.a. 

15. 

0.491 

0.230 

362 

16. 

XI • £L • 

n.a. 

310 

17. 

0.596 

0.275 

321 

18. 

0.6l4 

0.246 

228 

19. 

0.426 

0.336 

248 
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U.U Paxameters of the Second Type 

The geometrical parameters have all been nondimens ionalized through 
the use of the mean aerodynamic chord. The parameters are, therefore, 
defined as 


(U.3) 

b - c^/o 

(4.U) 

A. 

X^ = X^/c 

(4.5) 

A 

^2 ~ X 2 /C 

(4.6) 

X 3 = X^/c 

(4.7) 

/s 


^ ^ 

It should be noted that X^, X^, X^, and X^ represent distances from 
the craft's center of gravity to the location of the unsteady resultant 
forces acting on the wing and tail, as shown in Figure 2.1. They will 
be functions of spanwise position if the aircraft has a sweptback wing and 
tail and in such cases will be parameters of the third type. 

The characteristics of parameters of the second type are shown in 
Figure U.6. They have a constant value over the entire range of air- 
craft considered. 

The reduced frequency parameter for the wing is defined as 

0)C. 

(lie 

k = while that for the tail is similarly defined as k^ = . 

o o 

The ratio of the frequency parameters is, therefore, equal to the ratio 

of the mean aerodynamic chords shown in Figure U.6 to be 0.6. 
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GEOMETRICAL PARAMETERS 



Figure U,6 Computed Values for Geometrical Aircraft Parameters 
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4.3 Parameters of the Third Type 

Although assimed not to apply in the numerical analysis of this 
study, parameters of the third type occur due to change in a geometrical 
parameter vith spanwise location. 

The first case is that of a tapered wing and tail. The airfoil 
chord at any spanwise location may be expressed simply as the chord at 
the root multiplied by a coefficient a(y). The relation is therefore, 

(4.8) Cp(y) = a(y) 

The wing taper coefficient a(y) is best expressed in terms of 
another coefficient in the form, 

(4.9) a(y) = 1 - C^y 

where is expressible as, 

(4.10) 

Similarly for the tail, the following may be written. 

(4.11) C (y) = a (y) C 

^t ^ ""t 

(4.12) a^(y) = 1 - C^y 

lx - 1) 



The trend of the coefficients and are shown in Figure 4.7* 
Another contributing case to this third type of parameter is the swept 
wing case. In this case any distances from a fixed location on the 
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airfoil now become functions of y. As an example, the distances X^, X^, 

A A 

X^j and Xj^, which for a straight wing craft are invariant, now become 


(4.14) 

^ A 

X (y) = X 

j y tan A 1/4 


t 

c 

(4.15) 

X2(y) = X^ 

y tan A 1/4 


W 

c 

(4.16) 

/V A 1 

X (y) = X 

. y tan A 1/4 


W 

c 

(4.17) 

A A 1 

Xj^(y) = X^ 

y tan A 1/4 


The fundamental mode shape for elastic bending can be approximated 
by the expression 

(4.18) (|> 2 (y) = (l - cos 1^) 

and the fundamental torsional mode shape by 

(4.19) ~ if 
4.6 Exhibited Trends 

It is an extremely convenient and useful result that the various 
dimensionless parameters have such a simple form for conventional 
aircraft. This greatly generalizes the resulting equations of motion. 

In summary, the following characteristics are observed: Those 

parameters which depend on mass and mass distribution show a linear 
decreasing trend with increasing gross weight. Parameters which are 
geometrical without spanwise dependency are invariant over the entire 
weight range of craft. Geometrical parameters which are functions of 
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spanwise distance may be defined in terms of coefficients which are found 
to exhibit a linear decreasing trend with increasing gross weight. 
Parameters which are geometrical without spanwise dependence are invariant 
over the entire weight range of craft . 


k.l Bending to Torsion Frequency Ratio 

The equations of motion for elastic bending and torsion include the 
ratios of bending natural frequency to the forcing frequency and 

torsional natural frequency to the forcing frequency First, the 

value of the ratio of bending to torsion frequency will be estab- 

lished. 

This value may be obtained approximately by considering a cantilever 
wing. The expression for the elastic bending natural, frequency is 


( 1 +. 20 ) 



where c^^ is a coefficient whose value depends upon the mode considered. 
For the fundamental mode, its value is 3.5. 

The expression for the torsional natural frequency is 


( 1 +. 21 ) 



where again is dependent upon the torsional mode considered. 


The value of c for the fundamental mode is tt/2. 
a 


(U.22) 


The ratio of frequencies is, therefore, given by 


SII 


a 


GJmH 


The shear modulus G may be related to the Young's modulus E through 
the relation. 
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(it. 23) E = 2(1 + v)G 


where V is the Poisson ratio for the material, and is found to be about 
0.25 for aircraft materials. Therefore, 

(lt.2it) E = 2.5G 

The bending moment of inertia I is a geometrical quantity which may 
be expressed as : 

(it. 25 ) I = C^ct^ 

The other geometrical term in the frequency relation is the polar 
moment of inertia J which derives its primary contribution from bhe 
moment of inertia about an axis perpendicular to the chord line. There- 
fore, J is expressible as, 

(4.26) J = Cgtc^ 

The coefficient is approximately equal to coefficient for 
wing type sections. Since the ratio of the section moment of 

inertia may be expressed as, 

(4.27) (7) : 

J c 

The wings of subsonic jets have a thickness to chord ratio of 0.13» 
with few exceptions. 

The mass moment of inertia of a wing may be written as, 

(4.28) I = c mc^ 

o o 
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For a ving section the value of generally has a value of the 
order of 0.1. The ratio of c to Jl for subsonic jets has a value of O.UO. 
Therefore the following may be written: 

{h.29) i\) Z 0.1 

mZ 

Substituting equations (U.2U) through (4.29) into equation (4.22) 
yields : 

(4.30) — ; 1.13 (t/Z) 

0) 

a 

A reasonable approximation to the frequency ratio is, therefore: 

(4.31) — ::o.o6 
0) 

a 

Credibility is added to this value by an analysis performed by 
Houbolt and Anderson [l8] 19^8, who calculated the frequencies of the 
fundamental modes of a nonuniform beam. They found, using the Stodola 
method of frequency determination, that the bending frequency had the 
value in the fundamental mode of 2,h rad per sec. The fundamental 
torsional mode frequency is given as approximately 71*0 rad per second. 
This would indicate a value for the frequency ratio of 0.035, which is 
the same order of magnitude as is determined by this analysis. 

The effect of a variation in the frequency ratio 
established numerically. Computations were performed with a range 
of from 0.01 to 1.00 in the equations of motion. It was found 

that no shift in any resonant peak resulted. Slight variation in the 
meignitude did result, however, no consistent trend was verified and in all 
cases, the variation was less than ten percent of the magnitudes resulting 
from = 0 . 06 . 
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U.8 Relation of Natural Frequency to Aircraft Physical Parameters 


It would he advantageous to express the hending and torsional 
natxiral frequencies in terms of basic physical quantities and known 
aircraft parameters. Beginning with the relation for the bending 
natural frequency of a cantilever wing, the first mode is expressed, 

(4.32) u) 

h ^2 m 

After appropriate substitutions, the following relation is obtained: 


(U.33) 0^ = 1.T58 / (^) (|) If) (j|) 


The coefficients c^ arid Cg derive their values from the evaluation 
of the following expressions. The bending moment of inertia is expressed 


(4.34) I = 


c^ W At 


and the mass distribution is expressed as. 


(4.35) m = c^ pg tc 


Evaluation of coefficients in terms of those aircraft values 
which will be known for a given craft yields the expression for the 
coefficients as, 

O, I p c 


which may be written in terms of a coefficient previously introduced as. 
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(U.37) 


c’ ■ Wn 


The expression for frequency may also be stated in terms of air- 
craft parameters, values of which are given by Werner [l9l» Tbe 
expression for frequency is; 

/(§-) (S-)" (It-)' ( I) 


. t f t 


(U.38) u) = 9 X ^ 


c. ( 


c c 

V mx 


'a’ 'W ' 'ct^ 
b. s 


^6) 


Figures U.8 through U.IO give the Werner subsonic jet values . 

c c 3/2 . 

(U.39) c (-^5-) = 0.8 (10’^) 

a 


(U.Uo) ~ shown in Figure k.2) 

s 

(1..U) (4^) = (10-2) 


The value of (c/A) may be taken as O.UO for subsonic jet air- 
craft . 

The frequency for torsion may also be expressed in terms of 
several nondimensional parameters. Using the relation. 


(U.U2) 

O 


substitution of the following yields the required relation. 


(U.U3) J = 


c^’ W c^ 


s 


(U.UU) I = c’’ p c^t 
O 2 s 
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MAXIMUM TAKE-OFF WEIGHT, W x I0'^ LBS. 







The relation is 


(It. 1+5) w = 17.75 / U(^) 


a 




(£j) (f) (f) (^) 


A^ain the constants meiy be eval\iated from equations (4.U3) and (4,41+) . 

It may be seen therefore, that once the nondimens ional ratios of equations 
(It. 33) and (It. 1+5) are known, it is merely a matter of dictating a relative 
mass in order that the natural frequencies become known. 

The evaluation of the ratio yields the relation 

a p ct J 

(1..U6) 


In terms of the parameters as given by Werner, the relation for torsional 
frequency is. 


(It. 1+7) 0 )^ = cl' c ( 

Ot 3 C. 


3/2 / (V x3 /G . ,A .3 /Ax (B.) 


.^Using relations (1+.33) and (1+.1+5), the ratio of the frequencies 
may now he stated in the form. 


(It. 1+8) 


IlJh _ 

u " ^It*^ il SL 
a 


,1/2 


where A is the characteristic area of the aircraft. 

1/2 

For subsonic jet aircraft, both (t/Jl) and (A ' /I) are nearly 
constant over the entire range of craft. Therefore, it may be con- 
cluded that the ratio is a constant over the entire range. The 

value previously presented, equal to 0.06 is valid for all craft 

within the group. 
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5.0 NUMERICAL RESULTS AND CONCLUSIONS 


5.1 Power Spectral Techniques 

The equations formulated in Chapters 2.0 and 3.0 will yield the 
frequency response functions of the modes considered. For a stationary 
random process, power spectral techniques may he used to determine the 
power spectral density of output quantities once input spectral densities 
are established. Given a random input or output function, the power 
spectral density is defined as the Fourier transform of the auto- 
correlation function [ 25 ]. The power spectral density of the input may 
he related to the power density of the output hy the well known expression: 

( 5 . 1 ) 'F' (k) = |H(ik) 1^ $' (k) 

where ¥’(k) and $*(k) are the respective spectral densities. The response 
of the system depends on the input spectral density and secondly on the 
characteristics of the mechanical admittance squared. 

Mathematical expressions frequently used to represent power spectral 
density curves contain three parameters. The root mean square velocity 
measures the intensity of the turbulence. When the power spectral density 
is plotted against frequency, the r.m.s. has the effect of moving the 
curve upward or downward parallel to the ordinate. A second parameter 
determines the slope at the high frequency end of the curve. The third 
parameter, scale length, is a shape parameter defining the rate change 
of the slope in the low frequency region of the curve. 
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Relations for the input spectral density of atmospheric turbulence 
have been presented by numerous investigators. Von Karman [k] 19^0, 
presented an expression for the power spectral density of atmospheric 
turbulence based on a high frequency slope of-(5/3)- Crooks, Hoblit, 
and Prophet [32] 1967, have presented power spectral density relations 
appropriate to the case of high altitude clear air turb\ilence. Also, a 
study by Ashburn, Waco, and Mitchell [33] 1969, presents models of high 
altitude clear air turbulence. The scale length parameter distribution 
at high altitudes has been studied by Ashburn [3^] 1969- Most recently 
an expression for the input power spectral density has been presented by 
Houbolt [ 15 ], who presents aircraft structural design procedures based 
on this relation. The relation and its plot are shown in Figure 5*1* 

The purpose of this analysis is to present the characteristics 
of the mechanical admittance functions which generally depend on the 
relative mass parameters, Mach number, and the frequency at which the 
aircraft is being forced. All other aircraft parameters may either 
be related to one of these or are constant, as has been shown in 
Chapter 4.0. 

3.2 Mechanical Admittance Results 

For the class of subsonic transport jets considered in this analysis, 

the relative mass parameter has a range of two hundred (200) to six 

hundred ( 600 ). Equations (2.29), (2.30), (3.13), and (3.l4) maybe 

solved simultaneously to obtain the mechanical admittance values. If 

2 

equations (3.13) and ( 3 .lU) are divided by (k ), the ratio of bending 
to forcing reduced frequency (k^/k), and torsion to forcing reduced 
frequency (k^/k), occiir. In order to solve the equations, appropriate 
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values must te obtained for these ratios in terms of those quantities 
defined in Chapter U.O. 


The expression for torsionsil frequency was given previously as: 

(U.20) 


I I 

o 


The ratio of torsional frequency to forcing frequency may be 
expressed as follows: 


(5.2) 


a a /GJ 
k k I. 


where the coefficient c^ is given as 


(5.3) 


a 




The quantity a’ is the speed of sound at the flight altitude and M is the 
Mach number. 

In Chapter U.O, the expression relating the polar moment of 
inertia to the area moment of inertia about the chordwise axis was 
expressed as 

(1..27) (i) = (|)^ 

Using equations (U.25) and (U. 28 ), the ratio of the area moment 
of inertia to the mass moment of inertia becomes. 



o 


The quantity (J/I ), therefore, is written; 
o 

(5.5) (J/J ) ~ (~) 
o m 
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The nondimens ional wing mass distribution may be incorporated into 


equation (5-5) to yield; 


(5.6) 



t 

PmA 


Substitution of equation (5.6) into equation (5.2) yields for the re- 
quired ratio, 


(5.7) 



where c’’ is a function of altitude and Mach number, and may be written 

a 

as 

(5.8) =;■ = 0.067 /f 

Values for m, the nondimens ional wing mass distribution were given in 
Chapter k.O. 

The ratio of the bending reduced frequency parameter to the forcing 
reduced frequency may be stated: 

(5.9) (^) = f^) (^) 

a 

where determined to be approximately 0.06 in Chapter U.O. 

The results of manerical computations are tabulated in Appendix B 
for the cases of Mach number equal to 0.2, 0.5, and 0.7. Figures 5.2 
through 5.5 show the behavior of the mechanical admittances at a Mach 
number of O.J. 

Examination of Figures 5.2 through 5.5 show that secondary excita- 
tion occurs from modes other than the one being considered. In certain 
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instances, this secondary excitation is of nesurly comparable amplitude 
to that at the excitation frequency for that mode. This would indicate 
that the detrimental frequencies in a particular mode will be the char- 
acteristic frequency to which that mode responds plus possibly all 
frequencies characteristic to the other modes of the problem. 

The precision of the aerodynamic coefficients and other frequency 
dependent values could only be specified at intervals of reduced 
frequency equal to or greater than 0.01. Since both the plunging and 
pitching occur at relatively low frequencies, a highly detailed plot of 
each is impossible. The first resonant peak shown on Figures 5.2 through 
5.5 includes the contribution to the respective admittance quantities 
from both of these modes. If smaller intervel.s had been achievable, 
then the pitching and plunging modes could have been extracted separately. 
As shown in this scialysis, they appear together and no differentiation 
have been made between them. The critical values for both have been 
shown to be the same. 

5.3 Variation of Resonant Reduced Frequency with Mach Nximber and 

Relative Mass 

It is found that the resonant reduced frequencies remain nearly 
constant over the range of relative mass as shown in Figxires 5-6 through 
5.8. These values are taken from the niomerical results. A minor shift 
of the resonant reduced frequency to a slightly lower value does occur 
with increasing relative mass; however, the order of magnitude of the 
shift is less than * 0.005 . The effect of this shift is neglected and 
the frequency will be treated as invariant. 
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'igure 5*6 Resonant Reduced Frequency Variation with 
Relative Mass M = 0.2 
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RESONANT REDUCED FREQUENCIES (k) 



Figure 5*7 Resonant Reduced Frequency Variation with 
Relative Mass M = Q.5 
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Figure 5.8 Resonant Reduced Frequency Variation with 
Relative Mass M = 0.7 





The resonant reduced frequency is, however, dependent on the Mach 
number. This dependence shows a rapid growth of k as Mach number 
approaches zero. This would be expected from equation (5.3), which shows 
the frequency coefficient to be inversely proportional to Mach number. 

The variation of Mach number is illustrated in Figures 5*9 through 5.12. 
For Mach numbers of 0.2, 0.5, and 0.7, the resonant reduced frequencies 
have been established as shown in Table 5.1. The resonant frequencies 
for all Mach numbers are shown in Figures 5*13 through 5»l6 as derived 
from the numerical computation. 


Table 5.1 Resonant Reduced Frequencies 


Mach Number 

Plunging 

Pitching 

Bending 

Torsion 

0.2 

0.065 

0.065 

0.250 

3.290 

0.5 

0.025 

0.025 

0.098 

1.290 

0.7 

0.018 

0.018 

0.076 

1.000 


5.^ Variation of Admittance Magnitudes with Mach Number and Relative Mass 

Having considered the nat\ire of the resonant frequencies and their 
variation with Mach number and relative mass, the final question is how 
the magnitudes of the response are affected by variations in Mach number 
and relative mass. 

Figures 5.2 through 5.5 give the mechanical admittance values for 
the specific case of M = 0.7 and y = 200. It is found that the admittance 
functions for the plunging, pitching, and elastic bending modes decrease 
at y = 600 to approximately one third the value of the response at y = 200. 


63 




RESONANT REDUCED PLUNGING FREQUENCY Kr^ 



Figure 5*9 Variation of Resonant Plunging Reduced Frequency 
with Mach Number 
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REDUCED PITCHING FREQUENCY 



Figure 5.10 Variation of Resonant Pitching Reduced Frequency 
r with Mach Number 
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RESONANT RE 



Figure 5.11 Variation of Resonant Bending Reduced Frequency 
with Mach Number 
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RESONANT REDUCED TORSIONAL FREQUENCY 



MACH NUMBER. M 


Figure 5.12 Variation of Resonant Torsional Reduced 
Frequency with Mach Number 
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Figure 5*13 Variation of Resonant Circular Frequency 
(o)R^) with Relative Mass 
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Figure 5-15 Variation of Resonant Circular Frequency 

(co, ) with Relative Mass 
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Figure 5.16 Variation of Resonant Circular Frequency 
(a>^) with Relative Mass 
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These values of relative mass represent the extremes. Similarly, the 
torsional mode mechanical admittance decreases through this range to 
approximately one fourth of its value at y = 200. These decreases occur 
in a linear manner and. are depicted in Figures 5. IT through 5«32. 

The behavior of the mechanical admittance magnitudes with variations 
in Mach nxunber may be seen in Figures 5.33 through 5-36. Contrary to 
experience, the rigid body modes are found to be particularly sensitive 
to Mach nvunber variations. For this analysis, however, the controls 
have been assumed locked. This may account for the inconsistency 
between the niamericad results and experience. The rigid body modes 
increase in magnitude as the Mach number approaches zero due to the 
longer duration gust. The decreasing magnitude of the cycle frequency 
tends to reduce the amplitude of the elastic modes. This accounts for 
the contrast between the rigid and elastic modes. The magnitudes 
shown in the figures represent only the resonant magnitudes . 

5.5 Wavelength Critical Values for Resonant Response 

A study by Ashburn, Waco, and Melvin [3^], 1970, determined the 
wavelengths encountered in the high altitude clear air turbulence spectra. 
Root mean square gust velocities are established for the various values 
of wavelength and correlated to season of year, type of terrain, and 
altitude. It would be advantageous to determine those wavelength values 
which are critical to resonant response of an aircraft. 

The circular frequency of encounter for these wavelengths may be 
expressed as, 

2n u 

(5.10) to = 2. 

X 
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Figure 5.21 Variation of Admittance Magnitude H^(ik) with 
Relative Mass at Resonant Frequencies M = 0.2 
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.22 Variation of Admittance Magnitude H 2 (ik) with 
Relative Mass at Resonant Frequencies M =0.2 
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Figure Variation of Admittance Magnitude T^Cik) with 

Relative Mass at Resonant Frequencies M = 0.2 
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Figure 5*30 Variation of Admittance Magnitude H2 (ik) with 
Relative Mass at Resonant Frequencies M = 0.7 
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Figure 5*31 Variation of Admittance Magnitude T^(ik) with 
Relative Mass at Resonant Frequencies M = 0.7 
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Figure 5.33 Variation of Admittance Magnitude H^(ik) with 
Mach Number at Resonant Frequencies y = 200 
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Figure 5.3^+ Variation of Admittance Magniti 
Mach Number at Resonant Frequei 


90 



ADMITTANCE MAGNITUDE FOR T, (ik) 


1 



Figure 5-35 Variation of Admittance Magnitude T^(ik) with 
Mach Number at Resonant Frequencies jj = 200 
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The critical reduced forcing frequency may be e3q>r eased as 

(5.11) h = ^ 

X 

Each mode considered in the analysis will have a critical wavelength 
vhich will be dependent on the mean aerodynamic chord c , and the 
resonant reduced frequency. This value of wavelength is written: 

( 5 . 12 ) 

c 

The resonant reduced frequencies were found to be independent of 
the relative mass parameter. Nondimens ionali zing the resonant wave- 
length using the mean aerodynamic chord, the variation of the 
nondimens ional wavelength with Mach number is presented in Figures 5.37 
through 5.^0. 

The resonant wavelengths are plotted against relative mass of the 
craft, for Mach numbers equal to 0.2, 0.5* and 0.7» in Figures 5«^1 through 
5.U3. 

Conclusions 

The following conclusions may be extracted from the resvilts of this 
ansQysis. 

(1) The response of an aircraft to atmospheric turbulence is 
dependent on the parameters of relative mass, Mach number, and the reduced 
frequency at which it is being forced. 

(2) The ratio of mode reduced frequency parameter to forcing reduced 
frequency is dependent on Mach number and flight altitude. 
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Figure 5.k2 Resonant Turbulence Wavelength Variation with 
Relative Mass 
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Figure 5.^3 Resonant Turbulence Wavelength Variation 
with Relative Mass 


100 



( 3 ) Secondary excitation of inodes other than the node under con- 
sideration vill occur. The amplitude of secondary excitation can be 
comparable to the primary mode amplitude. 

(U) Resonant frequencies remain nearly constant vith variation in 
relative mass. These frequencies vary proportionally to the reciprocal 
of the Mach number at a given altitude. 

( ^ ) Tbe magnitudes of the rigid nodes increase vith decreasing Mach 
number while the elastic nodes decrease vith decreasing Mach Number.. 

( 6 ) The critical reduced frequencies are related to the reciprocal 
of the critical wavelengths of the turbulence spectra by the factor it 
times the mean aerodynamic chord. 

(7) The critical wavelengths of the turbulence spectra decrease in a 
logarithmically linear fashion with logarithmically increasing relative 
mass. 


5.7 Summery of Ma.jor Assumptions and Implications 

(1) The aircraft was assumed to be traveling in steady level flight 
at an altitude of thirty ’uiio'i isand feet. 

(2) Small disturbance approximations were made to linearize the 
response portion of the equations while the forcing terms were not 
linearized. This is a valid procedure as most disturbances do create 
response values of small magnitude. 

(3) The wing of the craft was assumed to be elastic. The remainder 
of the craft was assumed rigid. This is a reasonable engineering ap- 
proximation in considering only the pitching smd plunging rigid modes 
and the first bending and torsional modes. 
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(U) Linear trends were established for wing and aircraft mass 
dependent and geometrical parameters. These trends were only established 
for subsonic Jet transports and should be recomputed for other types of 
craft . 

(5) The wing was considered to be unswept and untapered. The 
existence of skew and variable cross section will alter frequency 
characteristics ajid forcing properties; however, the error introduced 
does not nullify the qualitative information deduced from the results. 

(6) wing elastic characteristics were established by treating the 
wing as a cantilever. Again, this introduces some error in the mode 
shapes assvuned, but is qualitatively a satisfactory engineering 
approximation . 

(7) The turbulence was assumed to be statistically stationary. 

Only vertical gusts were considered. 

5 . 8 Futxire Work 

An extension of this method could lead to increased accuracy 
and further insight into the effect of several factors on the response 
of a craft to random gusts. The following considerations are suggested: 

(1) The effects of spanwise variations of the unsteady aerodynamic 
loads on aircraft response. 

(2) The use of three-dimensional unsteady aerodynamic theories. 

(3) The incorporation of wing taper and sweepback. 

(4) The inclusion of additional bending and torsional elastic 
modes for the wing. 

(5) Consideration of other segments of the aircraft as elastic. 
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(6) The inclusion of chordwise gusts for analysis of the variations 
occurring in forward flight velocity. Lateral gusts might blLso he 
considered. 
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APPENDIX A 

A COMPUTER PROGRAM - 

numerical Evaluation of the admittahce values 
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A COMPUTER PROGRAM - 


NUMERICAL EVALUATION OF THE ADMITTANCE VALUES 


The numerical evaluation of the equations of motion may he 
accomplished through the use of the subroutine CLEQP which computes 
the solutions to a set of simultaneous linear equations with complex 
coefficients . 

The following program has been utilized in order to allow the 
variation of the relative mass parameter through its realistic range 
in these equations. 


IMPLICIT REAL* 8 (A.D-H, 0-Z) 

C0MPLEX*16 C(U0,4,4),C0M(U,U),B(1+),BC(4),BET 
REAL*8 MU(40) 

DIMENSION A(4,4,4) 

WRITE(6,105) 

105 F0RMAT( ) 

XMU=200 

DMU=20 

DO 3 IR0W=1,4 
DO 1 IC0L=1,4 

READ(5,100) (A(I,IR0W,IC0L),I=1,4) 

100 FORMAT ( 4 (F12.T,8X}) 

1 CONTINUE 
READ(5,102)B(IROW) 

102 FORMAT (2F20.T) 

3 CONTINUE 

DO 4 IMU=1,40 
MU(IMU)=XMU 
DO 2 IR0W=1,4 
DO 2 IC0L=1,4 

AX=A ( 1 , IROW , ICOL ) +A ( 2 , IROW , ICOL ) *XMU 
AY=A ( 3 , IROW , ICOL ) +A ( 4 , IROW , ICOL ) *XMU 
C UMU , IROW , ICOL ) =DCMPLX ( AX , AY ) 

2 CONTINUE 

4 XMU=XMU+DMU 
DO 5 IMU=1,40 
DO 6 IR0W=1,4 

BC ( IROW )=B( IROW) 

DO 6 IC0L=1,4 

6 COM (IROW, ICOL )=C(IMU, IROW, ICOL) 

CALL CLEQD(C0M,BC,4,1,4,4,BET) 
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WRITE(6,10U) MU(IMU) 
lOU F0RMAT('-',50X,‘M[J* •,F10.5/) 
WRITE(6,103) (BC(IR0W),IR0W=1,4) 

103 FORMAT{ ' Ansver= ’ , 5X ,D20 . 7 , 5X ,D20 . 7 , ’ *I ' ) 
5 COUTIHUE 
STOP 
END 
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APPEBDIX B 

MECHANICAL ADMITTANCE CCMffONENTS 


109 



MECHAUICAL ADMITTANCE COMPONENTS 


The following tables list the real and imaginary components of the 
mechanical admittance. These are given for a variety of significant 
reduced frequency values and cover the entire range of relative mass 
for subsonic jet transports. 
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Table B.l Mechanical Admittance Components 


Reduced Frequency ** 0.02 Mechanical Admittance Magnitude M = 0.2 


y 

H^(ik) 

Hgdk) 

T^(ik) 

T^dk) 

200 

.12716-02 

.17798-021 

O.20736-OU 

0.226 i 4-0U1 

-0.33522-04 

0.37468-041 

0.98736-06 

0.51647-071 

300 

O.8377U-O3 

0.13U02-021 

0.22817-OU 

0.17632-OUi 

-0.26149-04 

0.24729-041 

0.64337-06 

0.30849-071 

Uoo 

0.55061-03 

0.11019-021 

O.23035-OU 

0.13228-OUi 

. -0.21792-04 

0.17985-041 

0.44731-06 

0.21634-071 

500 

O.3U713-O3 

0.90U2U-03i 

0.21669-OU 

0.95173-051 

-0.20113-04 

'0.11990-041 

0.34768-06 

0.12976-071 

6oo 

0.22386-03 

0.76003-031 

0.19893-OU 

0.68172-051 

-0.18453-04 

0.93756-051' 

0.27369-06 

0.84316-081 

Reduced Frequency = 

O.oU Mechanical Admittance Magnitude M = 0.2 

y 

H^(ik) 

HgCik) 

T^(ik) 

TgCik) 

200 

-0.63612-02 

0.88993-021 

0.10513-03 

0.12631-031 

-0.10385-03 

-0.11858-031 

0.25642-04 

0.12637-051 

300 

-O.U1908-02 

0.67022-021 

O.II67U-O3 

-0.96932-OUi 

-0.79252-04 

-0.78543-041 

0.16849-04 

-0.73942-061 

Uoo 

-O.275UU-02 

0.55087-021 

0.11939-03 

-0.7369U-0U1 

-0.66955-04 

-0.56932-041 

0.11634-04 

-0.50732-061. 

500 

-0.17366-02 

O.U 5119-021 

Q. 10978-03 
-0.517UU-0U1 

-0.63001-04 

-0.38217-041 

0.93746-05 

-0.31769-061 

600 

« 

-0.11196-02 

0.38001-021 

0.98210-04 

-0.37422-041 

-0.57206-04 

-0.29554-041 

0.72632-05 

-0.3164 i- 06 i 
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TABLE B.l (Continued) 


Reduced Frequency = 0.06 Mechanical Admittance Magnitude M = 0.2 


y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.86702-01 

0.12135-OOi 

0.67559-03 

-0.74213-031 

-0.56950-03 

-0.63750-031 

0.39671-03 

-0.236i7-04i 

300 

-0. 57119-01 
0.91386-Oli 

0.75211-03 

-0.586l4-03i 

-0.44382-03 

-0.41583-031 

0.26379-03 

-0.i4736-04i 

iiOO 

-0.37542-01 

0.75131-Oli 

0.76843-03 

-0.45317-031 

-0.36855-03 

-0.30621-031 

0.18237-03 

-0.92632-051 

500 

-0.23670-01 

0.61528-Oli 

0.70966-03 

-0.30127-031 

0.34208-03 

-0.20483-031 

0.13614-03 

-0.69361-051 

6oo 

-0.15260-01 

0.51894-Oli 

0.63715-03 

-0.23144-031 

-0.15938-031 

-0.15938-031 

-0.42734-051 

-0.42734-051 

Reduced Frequency = 

0.065 Mechanical Admittance 

Magnitude M = 0.2 

y 

H^(ik) 


\(ik) 

Tgdk) 

200 

-0.40460-00 

0.56630-OOi 

0.30475-02 

0.33949-021 

-0.12067-01 

-0.13506-Oli 

0.13228-02 

-0.67336-041 

300 

-0.26655-00 

0.42648-001 

0.33742-02 

-0.26881-021 

-0.93614-02 

-0.89103-021 

0.86573-03 

-0.45239-041 

1+00 

-0.17531-00 

0.35060-001 

0.34020-02 

-0.21451-021 

-0.78074-02 

-0.64859-021 

0.59376-03 

-0.26933-041 

500 

-0.11046-00 

0.28711-OOi 

0.30893-02 

-0.14622-021 

0.72341-02 

-0.43214-021 

0.46321-03 

-0.6572-041 

600 

-0.71215-01 

0.24l8l-00i 

0.28316-02 

-0.10194-021 

0.66012-02 

-0.33757-021 

0.34270-03 

-0.11747-041 
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TABLE B.l (Continued) 


Reduced Frequency = 0.08 Mechanical Admittance Magnitude M = 0.2 


y H^(ik) H2(ik) T^(ik) T2(ik) 


200 

- 0 . 63580-02 

- 0 . 88992 - 02 i 

0 . 92664-03 

0 . 10649-021 

0.60355-03 

- 0 . 67531-031 

0.94765-04 

0.51637-051 

300 

- 0 . 41888-02 

- 0 . 67020-021 

0.10436-02 

0 . 82317-031 

0.47034-03 

- 0 . 44222 - 031 . 

0.65759-04 

0 . 30842-051 

Uoo 

-0.27549-02 

- 0 . 55097-021 

0.10756-02 

0.63187-031 

0.39195-03 

- 0 . 32649-031 

0.43662-04 

0.20696-051 

500 

-0.17357-02 

- 0 . 45119-021 

0.96037-03 

0.43632-031 

0.36018-03 

-0.21637-031 

0 . 34672-04 

0 . 13264-051 

600 

-0.11193-02 

- 0 . 380 l 4 - 02 i 

0.89371-03 

0.30026-031 

0.33165-03 

-0.16875-031 

0.26312-04 

0.91306-061 

Reduced Frequency = 

0 . 10 Mechanical Admittance 

Magnitude M = 0.2 

y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.26588-02 

0 . 37214-021 

0.72936-03 

0 . 82642-031 

0 . 41372-03 

-0.45376-031 

0 . 54639-04 

-0.26394-051 

300 

-0.17516-02 

-0.28026-021 

0.80562-03 

0 . 65147-031 

0.32659-031 

- 0 . 30744-031 

0.33947-04 

0.16543-051 

400 

-0.11521-02 

- 0 . 23039-021 

0 . 82043-03 

0.50477-031 

0.28552-03 

-0.23678-031 

0 . 24849-04 

0.10637-051 

500 

-0.72596-03 

-0.18867-021 

0.76265-03 

0.34717-031 

0 . 24479-03 

-0.16137-031 

0.17650-04 

0.62573-051 

6oo 

- 0 . 46797-03 

-0.15891-021 

0.69459-031 

0 . 25742-031 

0.22745-03 

-0.11632-031 

0 . 13705-04 

-0.43671-061 
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TABLE B . 1 ( Cont inued ) 


Reduced Frequency = 0.l6 Mechanical Admittance Magnitude M = 0.2 


M 

H^(ik) 

HgCik) 

T^(ik) 

Tgdk) 

200 

-0.15606-02 

0.21843-021 

-0.55673-03 

0.61575-031 

0.27564-03 

0.30611-031 

0.36214-04 

0.19735-051 

300 

-0.10281-02 

0.16450-021 

-0.6i410-03 

0.50731-031 

0.21745-03 

0.20967-031 

0.24379-04 

0.12598-051 

4oo 

-0.67638-03 

0.13523-021 

-0.62764-03 

0.37659-031 

0.18421-03 

0.15007-031 

0.16845-04 

0.81732-061 

500 

-0.42601-03 

0.11074-021 

-0.6o448-03 

0.25397-03 

0.16237-03 

0.97536-041 

0.12073-04 

0.50631-061 

6oo 

-0.27017-03 

0.93276-031 

-0.51243-03 

0.19551-031 

0.14035-031 

0.82714-041 

0.10078-04 

0.33649-061 

Reduced Frequency = 

0.22 Mechanical 

Admittance 

Magnitude M = 0 . 2 

y 

( ik ) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.23120-02 

0.32363-021 

-0.13574-02 

0.15173-021 

0.46592-03 

0.53171-031 

0.99331-03 ■ 
0.52997-041 

300 

-0.15231-02 

0.24373-021 

-0.14922-02 

0.12363-021 

0.35256-03 

0.38197-031 

0.65130-03 

0.31746-041 

400 

-0.10018-02 

0.20036-021 

-0.15273-02 

0.90455-031 

0.31778-03 

0.27335-031 

0.44384-03 

0.22374-041 

500 

-0.63124-03 

0.61409-021 

-0.14231-02 

0.63712-031 

0.27849-03 

0.17224-031 

0.35698-03 

0.12634-041 

6oo 

-0.40877-03 

0.13819-021 

-0.12659-02 

0.46128-031 

0.24639-03 

0.13004-031 

0.29173-03 

0.89312-051 
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TABLE B.l (Continued) 


Reduced Frequency = 0.25 Mechanical Admittance Msignitude M — 0.2 


y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.14450-00 

0.20225-00i 

-0.14362-01 

0.15873-011 

0.34785-03 

0.37914-031 

0.62364-01 

0.33419-021 

300 

-0.95203-01 

0.15231-00i 

-0.15936-01 

0.12964-011 

0.27145-03 

0.25871-031 

0.41183-01 

0.20445-021 

400 

-0.62617-01 

0.12521-OOi 

-0.16437-01 

0.95489-021 

0.23369-03 

0.19472-031 

0.27998-01 

0.13256-021 

500 

-0.39542-01 

0.10254-001 

-0.14842-01 

0.66765-021 

0.20843-03 

0.13791-031 

0.20172-01 

0.78463-031 

600 

-0.25436-01 

0.86361-Oli 

-0.13747-01 

0.48172-021 

0.18954-03 

0.96536-041 

0.16947-01 

0.52366-031 

Reduced Frequency = 0. 

50 Mechanical 

Admittance Magnitude M = 0.2 

y 

( ik ) 

H2(ik) 

T^(ik) 

T^dk) 

200 

-0.57802-03 

0.80936-031 

-0.24715-03 

0.27442-031 

-0.58342-04 

0.63772-041 

0.19973-04 

0.10647-051 

300 

-0.38080-03 

-0.60953-031 

-0.27236-03 

0.21547-031 

-0.46238-04 

0.42693-041 

0.13695-04 

0.63274-061 

400 

-0.25046-03 

-0.50107-031 

-0.27946-03 

0.15884-031 

-0.39528-04 

0.32147-04i 

0.93672-05 

0.43127-061 

500 

-0.15780-03 

-0.41035-031 

-0.28731-03 

0.11221-031 

-0.35714-041 

0.21935-041 

0.72679-05 

0.26334-061 

600 

-0.10173-03 

-0.34560-031 

-0.23514-03 

0.82241-031 

-0.32142-04 

0.16324-041 

0.61284-05 

0.17379-061 
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TABLE B.l (Continued) 


Reduced Frequency = 1.00 Mechanical Admittance Magnitude M = 0.2 



H^(ik) 

Hgdk) 

T^(ik) 

•H 

CM 

E^ 

200 

0.28913-03 

0.40450-03i 

-0.12367-03 

0.13649-031 

-0.23654-04 

0.27136-041 

0.14115-04 

0.73652-061 

300 

0.19048-03 

0.30463-03i 

-0.13647-03 

0.10748-031 

-0.l84l9-04 

0.15937-041 

0.96813-05 

0.4315^-061 

400 

0.12519-03 

0.25043-03i 

-0.13923-03 

0.83651-041 

-0.16431-04 

0.13613-041 

0.68317-05 

0.30105-061 

500 

0.78931-04 

0.20508-03i 

-0.12847-03 
0. 52761 -04i 

-0.14736-04 

0.92371-051 

0.49377-05 

0.18762-061 

600 

0.50899-04 

0.1T2T2-03i 

-0.11527-03 

0.4076i-04i 

-0.11844-04 

0.72779-051 

0.371^6-05 

0.12610-061 

Reduced Frequency = 2.00 Mechanical Admittance Magnitude M - 0.2 

y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.20231-03 

-0.28315-03i 

-0.87578-04 

0.98273-041 

-0.16598-04 

0.l8846-04i 

-0.12173-04 

-0.63179-061 

300 

-0.13328-03 

-0.21324-031 

-0.96277-04 

0.8l448-04i 

-0.12957-04 

0.12647-041 

0.83654-05 

-0.38762-061 

400 

-0.87603-04 

-0.17530-031 

-0.97863-04 

0.59314-041 

-0.10751-04 

0.94337-051 

0.54671-05 

-0.25992-061 

500 

-0 . 55234-04 
-0.14356-031 

-0.91443-04 

0.4o879-o4i 

-0.98563-05 

0.67336-051 

0.43612-05 

-0.15837-061 

6oo 

-0.35617-04 

-0.12091-031 

-0.83842-04 

0.28495-041 

-0.92634-05 

0.45316-051 

0.33674-05 

-0.10735-061 


U6 



T ABLE B.l’ (Continued) 


Reduced Frequency = 2.50 Mechanical Admittance M = 0.2 


y 

H^(ik) 

Hgdk) 

T^(ik) 

TgCik) 

200 

-0.3it 682-03 
-0.U85it0-03i 

-0.93613-Olt 

0.11375-031 

0.2l61t0-0U 

-0.22703-0Ui 

-O.12973-OU 

-0.6637it-06i 

300 

-O.228U9-O3 

-0.36555-031 

-0.10298-03 

0.88781t-0iti 

0.l6327-0it 

-0.15216-Olti 

-0.87364-05 

-0.1t0l75-06i 

Uoo 

-0.15028-03 

-0.30051-03i 

-O.IO6U1-O3 

0.659it3-0lti 

0.llt733-0lt 

-0.11536-OUi 

-0.59715-05 

-0.26798-061 

500 

-0.9it683-03 

-0.2U6l0-03i 

-0.9736it-0U 

0.1t7331-0lti 

O.12373-OU 

-0.7653it-05i 

-O.U3622-05 

-0.16327-061 

6oo 

-0.6ioU5-03 

-0.20T2T-03i 

-0.89137-Olt 

0.33752-Olti 

0.1llt70-0lt 

-0.55236-051 

-0.13254-05 

-0.11273-061 

Reduced Frequency = 3. 

QQ Mechanical Admittance Magnitude M = 0.2 

y 

H^(ik) 

HgCik) 

T^(ik) 

Tgdk) 

200 

-0.2lt276-02 

-0.33978-02i 

-0.18732-03 

0.20llt6-03i 

0.61731-Oit 

-0.70652-Olti 

-0.98736-03 

-0. 48124-04 1 

300 

-0.15993-02 

-0.25589-021 

-0.20613-03 

0.l61tlt6-03i 

0.1t7605-0lt 

-0.1t6937-0lti 

-0.64317-03 

-0.30120-041 

it 00 

-0.10519-02 

-0.21036-021 

-0.20937-03 

0.120l| 3-031 

0.1t0930-0lt 

-0.35101-Olti 

-0.44643-03 

-0.21316-04i 

500 

-0.66275-03 

-0.17227-021 

-0.19561-03 

0.91319-Oiti 

0.36128-Olt 

-0.2lt936-0lti 

-0.34138-03 

-0.12695-041 

600 

-0.U263it-03 
-0. lit 509-021 

-0.178iU-03 

0.62275-Oiti 

0.3261tl-0U 

-0.l8l3it-0iti 

-0.28312-03 

-0.83219-051 
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TABLE B.l (Continued) 


Reduced Frequency = 3.25 Mechanical Admittance Magnitude M = 0.2 


U H^Cik) HgCik) T^(ik) Tg(ik) 


200 

-0.26015-01 

-0.36U05-01i 

-0.10427-01 

0.11736-Oli 

0.63763-03 

-0.74126-031 

-0.14269-00 

-0.76413-021 

300 

-0.17139-01 

-0.274l7-01i 

-0.11932-01 

0.95178-021 

0.49134-03 

-0.49376-031 

-0.93947-01 

-0.48316-021 

Uoo 

-0.11272-01 

-0.22538-Oli 

-0.12174-01 

0.73459-021 

0.42607-03 

-0.37094-031 

-0.67390-01 

-0.31778-021 

500 

-0.71022-02 

-0.18547-Oli 

-0.10978-01 

0.48755-021 

0.36982-03 

-0.26139-031 

-0.48736-01 

-0.19210-021 

600 

-O.U 5017-02 
-0.15545-Oli 

-0.98746-02 

0.35473-021 

0.32337-03 

-0.18561-031 

-0.4o612-01 

-0.13147-021 

Reduced Frequency = 

3 . 50 Mechanical 

Admittance Magnitude M = 0.2 

U 

H^Cik) 

Hg(ik) 

T^(ik) 

Tgdk) 

200 

-0.23126-02 

-0.32360-021 

0.67336-03 

0.74238-031 

0.40844-04 

-0.45175-041 

-0.20347-02 

-0.94639-041 

300 

-0.15235-02 

-0.24373-021 

0.75168-03 

0.59233-031 

0.31643-04 

-0.30746-041 

-0.13657-02 

-0.59351-041 

4 00 

-0.10020-02 

-0.20034-021 

0.76814-03 

0.45117-031 

0.27956-04 

-0.22840-04i 

0.91783-03 

-0.381i4-04i 

500 

-0.63133-03 

-0.l6407-02i 

0.70842-03 
0 . 31247-031 

0.24263-04 

-0.15881-041 

-0.69317-03 

-0.24137-04i 

6oo 

-0.40707-03 

-0.13818-021 

0.64736-03 

0.22641-031 

0.21654-04 

-0.11367-041 

-0.63119-03 

=0.15204-04i 
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1!ABLE B.l (Continued) 


Reduced Frequency » 4.00 Mechsmiced Admittance Magnitude M = 0.2 


y 

H^Cik) 

HgCik) 

T^(ik) 

Tgdk) 

200 

-0.31796-03 

-0.44495-03i 

0.87694-04 

0.98l36-04i 

0.13568-04 

-0.15633-04i 

-0.19364-04 

-0.10148-051 

300 

-0.20947-03 

-0.33509-03i 

0.95768-04 

0.79139-04i 

0.98837-05 

-0.10364-04i 

-0.13412-04 

-0.6094i-06i 

400 

-0.13768-03 

-0.27547-031 

0.97367-04 

0.58971-04i 

0.89134-05 

-0.78337-051 

-0.98716-05 

-0.41368-061 

500 

-0.86800-04 

-0.22792-031 

0.91473-04 

0.4l788-04i 

0.82155-05 

’-0.53214-051 

-0.68445-05 

-0.26198-061 

600 

-0.55963-04 

-0.18999-031 

0.83212-04 

0.29l47-04i 

0.68493-05 

-0.41369-051 

-0.60310-05 

-0.16217-061 

Reduced Frequency = 10.00 Mechanical Admittance Magnitude M = 0.2 

y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.11560-03 

-0.l6l80-03i 

0.41378-04 

0.45H3-04i 

0.68314-05 

-0.75634-051 

-0.64397-05 

-0.33649-061 

300 

-0.76165-04 

-0.12185-031 

0.45367-04 

0.35942-04i 

0.48813-05 

-0.50158-051 

-0.43164-05 

-0.20938-061 

400 

-0.50099-04 

-0.10017-031 

0.46229-04 

0.27l48-04i 

0.45137-05 

-0.38426-051 

-0.27931-05 

-0.13597-061 

500 

-0.31565-04 

-0.82039-041 

0.43237-04 

0.l8762-04i 

0.42368-05 

-0.25993-051 

-0.23102-05 

-0.84910-071 

600 

-0.20014-04 

-0.69093-04i 

0.40276-04 

0.13647-04i 

0.35362-05 

-0.19338-051 

-0.17604-05 

-0.64372-071 
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TABLE B . 1 ( Cont inued ) 


Reduced Frequency = 0.02 Mechanical Admittance Magnitude M = 0.5 


y 

H^(ik) 

HgCik) 

T^(ik) 

Tgdk) 

200 

0.46240-01 

0.64724-Oli 

0.42617-03 

0.38126-031 

-0.38573-05 

0.22647-051 

0.98437-04 

0.12642-041 

300 

0.30472-01 

0.48737-Oli 

0.25195-03 

0.50623-031 

-0.25679-05 

0.13264-051 

0.11869-03 

0.31645-04i 

1+00 

0.20021-01 

0.40064-01i 

0.15768-03 

0.50162-031 

-0.22364-05 

0.80126-061 

0.12439-03 

0.32243-041 

500 

0.12623-01 

0.32815-Oli 

0.11472-03 

0.49294-031 

-0.19714-05 

0.55273-061 

0.12314-03 

0.28136-041 

6oo 

0.81382-02 

0.27637-Oli 

0.92237-04 

0.48410-031 

-0.17264-05 

0.32419-06i 

0.11629-03 

0.25364-041 

Reduced Frequency = 

0.025 Mechanical 

Admittance Magnitude M = 0 . 5 

y 

H^(ik) 

H2(ik) 

T^(ik) 

Tgdk) 

200 

-0.18207-00 

0.25483-00i 

0.76437^02 

0.64169-021 

-0.73641-02 

-0.41665-021 

0.34852-02 

0.44219-031 

300 

-0.11998-00 

0.19189-OOi 

0.45331-02 

-0.85720-021 

-0.50175-02 

-0.22943-021 

0.41337-02 

-0.11039-021 

400 

-0.78836-01 

0.15774-OOi 

0.31449-02 

-0.84935-021 

-0.44320-02 

-0.14362-021 

0.43621-02 

-0.11845-021 

500 

-0.49705-01 

0.12920-00i 

0.20324-02 

-0.82736-021 

-0.37103-02 

-0.10638-021 

0.42932-02 

-0.98361-031 

6oo 

-0.32044-01 

0.1088l-00i 

0.16362-02 

-0.80213-021 

-0.31639-02 

-0.65741-021 

0.4i436-02 

-0.89445-031 
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TABLE B.l (Continued) 


Reduced Frequency = O.OU Mechanical Admittance Magnitude M = 0.5 


y 

H^(ik) 

H^Uk) 

T^(ik) 

Tgdk) 

200 

-0.69260-03 

0 . 97080 - 03 i 

0.3691+3-03 

0 . 32 ll+ 9 - 03 i 

-0.19873-01+ 

-0.11570-Ol+i 

0.12633-03 

0 . 14362-041 

300 

-0.1+5708-03 

0 . 73101 - 03 i 

0.22358-03 

-0.1+3627-031 

-O.I3572-OI+ 

-0.62736-051 

0 .i 4035-03 

- 0 . 35921-041 

Uoo 

-0.30033-03 

0 . 60093 - 03 i 

O.II+935-O3 

- 0 . 1 + 2886-031 

-O.I2OII+-OI+ 

-0.38975-051 

0.15367-03 

- 9 . 36575-041 

500 

-0.18935-03 

0 . 1 + 9220-031 

0 . 1061 + 3-03 

- 0 . 1 + 101 + 5-031 

- 0 . 961 + 27-05 

- 0 . 28413-051 

0 . 15144-03 

- 0 . 3214 o -041 

6oo 

-0.12207-03 

0 . 1 + 11 + 53-031 

O.6237I+-OI+ 

- 0 . 1 + 0310-031 

-0.87510-05 

-0.20078-051 

0 . 14132-03 

- 0 . 29436-041 

Reduced Frequency = 0.06 

Mechanical 

Admittance Magnitude M = 0.5 

y 

( ik ) 

Hgdk) 

T^(ik) 

T2(ik) 


200 

- 0 . 46240-03 

0 . 64720-031 

0.28374-03 

-0.23631-031 

-0.87253-05 

-0.49362-051 

0.10132-03 

- 0 . 12289-041 

300 

- 0 . 30472-03 

0 . 48734-031 

0.17295-03 

- 0 . 31268-031 

-0.60139-05 

-0.30227-051 

0.11963-03 

- 0 . 30417-041 

400 

-0.20022-03 

0 . 40062-031 

0.11296-03 

- 0 . 30847-031 

- 0 . 52461-05 

-0.17360-051 

0 . 12483-03 

- 0 . 31249-041 

500 

-0.12623-03 
0 . 32813-031 

0.83259-04 

-0.29654-031 

0 . 44219-05 

- 0 . 12642-051 

0.12361-03 

- 0 . 28436-041 

600 

-0.81382-04 

0.27635-031 

0.62368-04 

-0.28232-031 

0.39265-05 

-0.89921-061 

0 . 11845-03 

- 0 . 25147-041 
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TABLE B.l (Continued) 


Reduced Frequency = 0.08 Mechanical Admittance Magnitude M = 0.5 


y 

H^(ik) 

H2(ik) 

T^(ik) 

T^dk) 

200 

-0.80920-03 

-0.11326-02i 

0.33678-03 

0.29620-031 

0.98425-05 

-0.62194-051 

0.14236-03 

0.l6462-04i 

300 

-0.53313-03 

-0.85285-03i 

0 . 20147-03 
0.40219-03i 

0.66553-05 

-0.34632-051 

0.16843-03 

0.40687-04i 

4oo 

-0.35030-03 

-0.70107-03i 

0.12645-03 

0.39642-031 

0.59142-05 

-0.21258-051 

0.18269-03 

0.4l673-04i 

500 

-0.22091-03 

-0.57’+22-03i 

0.83629-04 

0.38126-031 

0.49137-05 

-0.15654-051 

0.17891-03 

0 . 37780-041 

600 

-O.1I+2U2-O3 

-0.1t8362-03i 

0.62241-04 

0.36764-031 

0.43177-05 

-0.11267-051 

0.40571-03 

0.33438-041 

Reduced Frequency = 0.10 

Mechanical 

Admittance Magnitude M = 0.5 

y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.60690-01 

0.81t9l*5-01i 

-0.37912-01 

0.32551-Oli 

0.22014-02 

-0.12643-021 

0.60377-01 

-0.72639-021 

300 

-0.39988-01 

-0.6396U-01i 

0.24622-01 

0.42659-Oli 

0.14731-02 

-0.739^+8-031 

0.73614-01 

0.18445-011 

Uoo 

-0.26279-01 

-0.52581-Oli 

0.14362-01 

0.41238-011 

0.13240-02 

-0.43655-031 

0.75267-01 

0.18927-Oli 

50(5 

-0.16568-01 

-0.1»3067-01i 

0.96357-02 

0.39645-011 

0.11462-02 

-0.31576-031 

0.74632-01 

0.159^+5-011 

6oo 

-0.10681-01 

-0.36271-Oli 

0.71430-02 

0.37845-011 

0.99873-03 

-0.21437-031 

0.73412-01 

-0.14501-011 
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TABLE B.l (Continued) 


Reduced Frequency = 0.12 Mechanical Admittance Magnitude M = 0.5 


y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.68377-03 

0.961»73-03i 

0.39213-03 

0.34621-031 

0.22147-04 

-0.12563-041 

0.59334-03 

-0.71462-041 

300 

-O.U5060-03 

-0.726Ult-03i 

0.24651-03 

0.45972-031 

0.i4639-04 

-0.70450-051 

0.72365-03 

0.17955-031 

Uoo 

-O.2956U-O3 

-0.59715-03i 

0.16672-03 

0.44867-031 

0.13562-04 

-0.41651-051 

0.75121-03 

0.18366-031 

500 

-0.18667-03 

-0.1t8912-03i 

0.10236-03 

0.43659-031 

0.11072-04 

-0.30651-051 

0.74o41-03 

0.15614-031 

600 

-0.12034-03 

-0.Ull94-03i 

0.86517-04 

0.41360-031 

0.99362-05 

-0.22641-051 

0.72109-03 

-0.14836-031 

Reduced Frequency = 0.l6 

Mechanical 

Admittance Magnitude M = 0.5 

y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.26010-03 

0.36405-03i 

-0.21635-03 

0.20146-031 

0.52639-05 

0.30120-051 

0.85133-04 

0.10372-041 

300 

-0.17141-03 

0.274l3-03i 

-0.13218-03 

0.28462-031 

0.41362-05 

0.17369-051 

0.10237-03 
0. 26130-04 1 

1*00 

-0.11262-03 

0.22535-03i 

-0.86342-04 

0.27659-031 

0.32490-05 

0.11225-051 

0.10784-03 

0.26845-041 

500 

-0.71007-04 

0.l8547-03i 

-0.62569-04 

0.26836-031 

0.26349-05 

0.75631-061 

0.10625-03 

0.22643-041 

6 oo 

-0.45778-04 

0.15545-031 

-0.45661-04 

0.25167-031 

0.22614-05 

0.52337-061 

0.10361-03 

0.21043-041 
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TABLE B.l (Continued) 


Reduced Frequency = 0.50 Mechanical Admittance Magnitude M = 0.5 


y 

H^(ik) 

Hg(ik) 

T^(ik) 

Tgdk) 

200 

-0.12138-03 

0.l6989-03i 

-0.98463-04 

0.87912-04i 

-0.10847-05 

0.67312-061 

0.51446-04 

0.62639-051 

300 

-0.79989-oU 

-0.12793-03i 

-0.57324-04 

0.11436-031 

-0.69147-06 

0.41362-061 

0.60395-04 

0.15637-041 

1+00 

-0. 52557-01+ 
-0.105l6-03i 

-0. 30463-04 
0.10592-031 

-0.60149-06 

0.24632-061 

0.63774-04 

0.15946-041 

500 

-0.33137-04 

-0.86l34-04i 

-0.23659-04 

0.98746-04i 

-0.51493-06 

0.17241-061 

0.63142-04 

0.13841-041 

600 

-0.21363-04 

-0.72543-04i 

-0.17931-04 

0.97240-04i 

-0.46107-06 

0.11579-061 

0.60226-04 

0.12667-041 

Reduced Frequency = 0. 

70 Mechanical 

Admittance Magnitude M = 0.5 

y 

H^(ik) 

Hgdk) 

T^^dk) 

Tgdk) 

200 

-0.86704-04 

0.12135-031 

-0.91778-04 

0.80937-04i 

0.10639-05 

-0.64172-061 

0 . 55372-04 

0.71073-051 

300 

-0.57138-04 

-0.91376-041 

-0.55935-04 

0.10639-031 

0.71632-06 

-0.35216-061 

0.66681-04 

0.17743-041 

Uoo 

-0.37543-04 

-0.75115-04i 

-0.35662-04 

0.10132-031 

0.62147-06 

-0.23798-061 

0.68938-04 

0.18013-041 

500 

-0.23670-04 

-0.6l524-04i 

-0.25184-04 

0.94713-031 

0 . 53692-06 
-0.16251-061 

0.67845-04 

0.15936-041 

600 

-0.15260-04 

-0.5l8l6-04i 

-0.19627-04 

0.89217-031 

0.45620-06 

-0.10639-061 

0.65840-04 

0.14324-041 
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TABLE B.l (Continued) 




Reduced Frequency = 1.00 Mechanical Admittance Magnitude M = 0.5 


y 


H2(ik) 

T^(ik) 

Tgdk) 

200 

O.I73U2-O3 

0.2U270-03i 

-0.11758-03 

0.98U36-0hi 

-0.15635-05 

0.85637-061 

0.101*1*2-03 

0.12552-01*1 

300 

0.1ll^28-03 

0.l8275-03i 

-0.68h39-Olt 

0.130U6-031 

-0.1121*0-05 

0.1*6731-061 

0.12837-03 

0.3081*7-01*1 

ItOO 

O.75O9I-OU 

0.15231-031 

-0.1»2108-04 

0.129llt-03i 

-0.90362-06 

0.281*76-061 

0.12681*-03 

0.31256-01*1 

500 

O.U73UU-OU 

0.12305-031 

-0.30l»52-0h 

0.11378-031 

-0.7926U-06 

0.22361-061 

0.121*79-03 

0.281*13-01*1 

600 

O.3052I-OU 

0.10363-031 

-O.2U337-OU 

0.10225-031 

-0.621*37-06 

0.15235-061 

0.12126-03 

0.25361-01*1 

Reduced Frequency =1.20 Mechanical 

Admittance Magnitude M = 0.5 

V 


HgCik) 

T^(ik) 

Tg(lk) 

200 

-0.85891-03 

-0.12556-021 

-0.1181*9-02 

0.1036U-021 

-0.701*36-05 

0.1*1698-051 

O.7968U-OI 

-0.10362-011 

300 

-0.56602-03 

-0.9i^5U6-03i 

-0.711*75-03 

0.1321*9-021 

_O.J*9175-05 

0.2339l*-05i 

0.96337-01 

-0.25638-011 

hOO 

-0.37195-03 

-0.77722-031 

-0.1*51*66-03 

0.1261*7-021 

-0.1*2357-05 

0.11*162-051 

0.1019l*-00 

-0.26238-011 

500 

-O.23U8-O3 

-0.63658-031 

-0.31*102-03 

0.11236-021 

-0.36113-05 

0.10971-051 

0.99235-01 

-0.22682-011 

600 

-0.15117-03 

-0.5361U-031 

-0.2361*1-03 

0.101*75-021 

-0.30759-05 

0.661*88-061 

0.95838-01 

-0.2103l*-015 




TABLE B.l (Continued) 


Reduced Frequency =1.30 Mechanical Admittance Magnitude M = 0.5 


u 

H^(ik) 

Hgdk) 

T^(ik) 

T2(ik) 

200 

-0.11097-01 

-0.15209-01i 

-0.25163-01 

0.2l65l*-01i 

-0.37521-01* 

0.21792-01*1 

0.45163-00 

-0.51320-011 

300 

-0.73129-02 

-0.1lL52-01i 

-0.15362-01 

0.2791*6-011 

O.25273-OI* 

0.12l61*-0l*i 

0 . 54265-00 
-0.12637-OOi 

400 

-0.1*8050-02 

-0.9l*ll*3-02i 

-0.82164-01 

0.26751-Oli 

-0. 21976-01* 
O.78I* 34-051 

0.56377-00 

-0.13149-001 

500 

-0.30295-02 

-0.77081+-02i 

-0.66370-01 
0.2 5228-01 i 

-0.19123-04 

0.53162-051 

0.55438-00 

-0.11617-OOi 

6oo 

-0.19531-02 

-0.6U9l*2-02i 

-0.1*3217-01 

0.2l*103-01i 

-0.17361-04 

0.36413-051 

0.54175-00 

-0.10436-001 

Reduced Frequency = l.UO 

Mechanical 

Admittance Magnitude M = 0.5 

y 

H^(ik) 

H2(ik) 

T^(ik) 

T2(ik) . 

200 

-0.11502-02 

-0.l6633-02i 

0.15632-02 

0.1361*3-021 

0.71352-05 

-0.40371-051 

-0.79972-02 

-0.10362-021 

300 

-0.75798-03 

-0.12525-02i 

0.93625-03 

0.17855-021 

0.47154-05 

-0.23472-051 

-0.10365-01 

-0.25372-021 

ItOO 

-0.1*9804-03 

-0.10296-021 

0.52631-03 

0.171*53-021 

0.41934-05 

-0.14361-051 

-0.10938-01 

-0.26297-021 

500 

-0.311*00-03 

-0.81*329-031 

0.1*1639-03 

0.16652-021 

0.35253-05 

-0.10362-051 

-0.10744-01 

-0.22648-021 

6oo 

-0.2021*3-03 

-0.71023-03i 

0.32658-03 

0.15297-021 

0.30567-05 

-0.65343-061 

-0.10298-011 

-0.20561-021 
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TABLE B.l (Continued) 


Reduced Frequency = 2.00 Mechanical Admittance Magnitude M = 0.5 


y 


HgUk) 

Tj^(ik) 

Tgdk) 

200 

- 0 . 68207 -Olt 

- 0 . 99917 -OUi 

0 . 7639 **- 0 U 

0.65115-OlH 

0.78527-06 

- 0 . 1 t 6 l 04 - 06 i 

-O.UU839-OI* 

-0.55178-051 

300 

- 0 . 1 tU 948 “ 0 l* 

- 0 . 75230 - 0 Ui 

O.U5693-OU 

0.87212-oUi 

0.51*632-06 

-0.27lU9-06i 

-0.53620-01* 

- 0 . 1362 l*- 0 l*i 

400 

- 0 . 29534 -OU 

- 0 . 6 l 9 Ul- 0 Ui 

0 . 30975 - 0 k 

0 . 85932 - 0 l»i 

O.U7256-06 

-0.16705-061 

-0.56217-01* 

- 0 . 13918 -OIH 

500 

-0.18621-0U 

- 0 . 50658 - 0 Ui 

0 . 200 U 3 - 0 l» 

0 . 83269 - 0 l*i 

0.39620-06 

-0.12127-061 

-O.55876-OU 

- 0 . 1 l 87 l*- 0 l*i 

6oo 

- 0 .l 200 lt- 0 U 

-0.U266U-0lti 

O.16228-OI1 

0.8l3U2-0Ui 

0.31*053-06 

- 0 . 781 * 85-071 

-O.5332I-OI* 

-O.llOlU-OUi 

Reduced 

Frequency = 2.50 

Mechanical 

Admittance Magnitude M = 0.5 

y 

H^(ik) 

HgCik) 

T^(ik) 

Tgdk) 

200 

- 0 . 52023 -OU 

-0.70796-OUi 

0 . 60327 - 0 l^ 

0 . 5 U 639 - 0 l»i 

0.62819-06 

- 0 . 36 l 5 l*- 06 i 

-0.1*0321-01* 

- 0 . 50268-051 

300 

- 0 . 3 U 283 - 0 k 

- 0 . 533 O 9 -OUi 

O.382UI-0U 

0.70229-OUi 

0.1*3673-06 

-0.2051*8-061 

-O.U8136-OU 

- 0 . 1263 l*- 0 l*i 

Uoo 

-O.22526-0U 

-0.43823-OUi 

0.23669-Ok 

0.69l3U-0ki 

0.37285-06 

-0.13781-061 

-O.U9638-OU 

’- 0 . 12879 - 01*1 

500 

- 0 .llt 202 - 0 lt 

- 0 . 3589 U- 0 Ui 

0.17276-Olt 

0.6783U-0Ui 

0.31628-06 

-0.98716-071 

-0.1*8936-01* 

-0.11657-01*1 

600 

-0.91560-05 

- 0 . 30230 - 0 ki 

O.13173-0U 

0.66215-oUi 

0 . 2733 l *-06 

-0.6201*8-071 

-0.1*791+8-01* 

-0.10062-01*1 
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TABLE B . 1 (Continued) 


Reduced Frequency = 0.02 Mechanical Admittance Magnitude M = 0.7 


y 

H^(ik) 

H2(ik) 

T^(ik) 

T2(ik) . 

200 

0.79856-02 

0.11393-01i 

0.39492-02 

0.31030-021 

-0.14691-02 

0.16518-021 

0.55262-03 

0.11958-031 

300 

0.78207-02 

0.12946-01i 

0.23731-02 

0.17433-021 

-0.11470-02 

0.10955-021 

0.35081-03 

0.69724-041 

4oo 

0.72144-02 

0.l4383-01i 

0.15113-02 

0.10741-021 

-0.99369-03 

0.77842-031 

0.24693-03 

0.44783-041 

500 

0.62281-02 

0.15719-Oli 

0.98072-03 

0.64177-031 

-0.90230-03 

0.55989-031 

0.18376-03 

0.27978-041 

6oo 

0.48859-02 

0.l6886-01i 

0.63819-03 

0.31470-031 

-0.83550-03 

0.39081-031 

0.l4l85-03 

0.i4532-04i 

Reduced Frequency = 0*0U 

Mechanical 

Admittance Magnitude M = 0.7 

y 


Hgdk) 

T^(ik) 

T^dk) 

200 

-0.17358-02 

0.4l435-02i 

0.72182-02 

0.15006-021 

-0.12636-03 

-0.12202-031 

0.75615-03 

0.84050-04i 

300 

-0.19740-02 

0.43312-021 

0 . 54771-02 
-0.21233-031 

-0.10068-03 

-0.11066-031 

0.52024-03 

-0.59133-051 

4 00 

-0,24861-02 

0.4l54l-02i 

0.44567-02 

-0.94187-031 

-0.72640-04 

-0.12230-031 

0.398i4-03 

-0.45829-041 

500 

-0.29947-02 

0.36643-021 

0.38788-02 

-0.12356-021 

-0.4i473-04 

-0.13096-031 

0.32910-03 

-0.64i33-04i 

600 

-0.33110-02 

0.29925-02 

0.35198-02 

-0.12894-021 

-0.11717-04 

-0.13011-031 

0.28649-03 

-0.698i6-04i 
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TABLE B.l (Continued) 


Reduced Frequency = 0.06 Mechanical Admittance Magnitude M = 0.7 


y 

H^(ik) 

Hg(ik) 

Tj^dk) 

Tgdk) 

200 

-0.U2285-02 

0.68309-021 

0.71873-02 

-0.18318-021 

-O.3520U-O3 

-0.28100-031 

O.3U122-O3 

-0.U1386-0U1 

300 

-O.57U3U-O2 

0.U7592-02i 

0.98000-02 

-0.18067-021 

-0.12573-03 

-0.28U25-031 

O.U0299-O3 

-0.10222-031 

Uoo 

-O.5696U-O2 

0.30550-021 

0.10195-01 

-0.1UU5U-021 

-0.18915-OU 

-0.22210-031 

O.U3523-O3 

-0.10UU5-031 

500 

-O.51U0U-O2 

0.199UU-021 

0.95330-02 

-0.121U3-021 

-0.25673-OU 

-0.16561-031 

O.U2802-O3 

-0.9UU59-0U1 

600 

-O.U528U-O2 

0.13606-021 

0.86301-02 

-0.10887-021 

O.U3258-OU 

-0.12U78-031 

0.U0311-03 

-0.8U556-0U1 

Reduced Frequency = 0. 

.08 Mechanical Admittance Magnitude M = 0.7 

y 

H^(ik) 

Hgdk) 

Tj^dk) 

Tgdk) 

200 

-O.930U9-O2 

-0.3U686-0Ui 

0.16153-01 
0 .136^5-011 

0.11725-OU 

-0.U63U1-031 

0.10992-02 

0.52227-031 

300 

-0.51556-02 

-0.2676U-021 

0.93289-02 

0.17871-Oli 

0.12633-03 

-0.18U58-031 

O.7163U-O3 

0»75880-031 

Uoo 

-0.30796-02 

-0.2U67U-02i 

0.59206-02 

0.17U59-011 

0.10067-03 

-0.8U026-0U1 

O.U966O-03 

0.76480-031 

500 

-0.20999-02 

-0.20600-02i 

* O.U29U2-O2 
0.16773-Oli 

0.76617-OU 

-0.U53U2-0U1 

0.37848-03 

0.74290-031 

600 

-0.15631-02 

-0.17295-021 

O.338U8-O2 

0.16221-Oli 

0.60289-OU 

-O.27U03-OU 

0.30657-03 

0.72185-031 
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TABLE B.l (Continued) 


Reduced Frequency = 0.10 Mechanical Admittance Magnitude M = 0.7 


y 

H^(ik) 

HgUk) 

T^(ik) 

Tg(ik) 

200 

-0.23296-02 

0.108T9-03i 

0.73786-02 

0.80919-021 

O.I380U-O3 

-0.76615-04 i 

O.13I+77-O2 

-0.71010-01+1 

300 

-0.20256-02 

-0.75l+52-0l+i 

0. 811+58-02 
0.61+850-021 

0.10211-03 

-O.l+2908-Ol+i 

O.III+5O-O2 

0.1+3l+2l+-0l+i 

Uoo 

-0.17 81+2-02 * 
-0.86613-Ol+i 

O.81I+7I-O2 

0.1+7728-021 

0.809^0-01+ 

-0.26903-01+1 

0.10089-02 

0.39569-Ol+i 

500 

-0.15787-02 

-0.7081+8-0l+i 

O.7690U-O2 

0.33812-021 

0.67115-01+ 

-0.18091-Ol+i 

O.893I+7-O3 

0.11318-Ol+i 

6oo 

-0.11+036-02 

-0.39685-Ol+i 

0.70612-02 

0.23653-021 

0.57U01+-0U 

-0.12800-01+i 

0.79328-03 

-O.ll+81+2-Ol+i 

Reduced Frequency = 0.12 

Mechanical Admittance Magnitude M = 0.7 

y 

H^(ik) 

Hgdk) 

T^(ik) 

TgUk) 

200 

-0.18772-02 

0.1066l-02i 

-0.90329-03 

0.52585-021 

0.81006-0I+ 

-0.32235-Ol+i 

0.25395-03 

0.891+1+1+-01+1 

300 

-0.15513-02 
' 0.69970-031 

-0.1871+9-02 

0.1+2781+-021 

O.70328-0I+ 
-0. 11+81+ 1+-01+1 

O.IO8I+I+-O3 

0.68126-Ol+i 

1+00 

-0.13303-02 

0.1+8858-031 

-0.22072-02 

0.3l+08i+-02i 

0.60I+20-0I+ 

-0.71915-051 

O.I+725O-OI+ 

O.l+l+167-Ol+i 

500 

-0.11622-02 

0.35231-031 

-0.22639-02 

0.27391-021 

O.52525-OI+ 

-0.301+75-051 

O.I86IO-OI+ 

0.2631+7-Ol+i 

6oo 

-O.IO28I+-O2 

0.26019-031 

-0.22038-02 

0.2237U-021 

O.I+62I+5-OI+ 

-0.59567-061 

O.I+312I+-O5 

O.ll+059-Ol+i 
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HA£L£ B.l (Continued) 


Reduced Frequency = 0.l6 Mechanical Admittance Magnitude M = 0.7 


u 

H^(ik) 

Hgdk) 

T^dk) 

Tgdk) 

200 

-0.10851-02 

0.67026-031 

-0.65192-03 

0.28731-021 

0.45289-04 

0.10332-041 

0.29555-03 

0.15432-041 

300 

-0.90257-03 

0.1+2788-031 

-0.98139-03 

0.20709-021 

0.34265-04 

0.10339-041 

0.18479-03 

0.90165-051 

UOO 

-0.76263-03 

0.28982-031 

-0.10122-02 

0.151+61-021 

0.27459-04 

0.96001-051 

0.12954-03 

0.58165-051 

500 

-06.51+67-03 

0.2061+3-031 

-0.91+396-03 

0.12029-021 

0.22837-04 

0.87955-051 

0.96650-04 

0.38670-051 

6oo 

-0.57068-03 

0.15327-031 

-0.881+25-03 

0.96958-031 

0.19504 -o 4 

0.80482-051 

0.74031-04 

0.25710-051 

Reduced Frequency = 0.20 

Mechanical 

Admittance Magnitude M = 0.7 


H^(ik) 

Hgdk) 

T^dk) 

Tgdk) 

200 

-0.70618-03 

0.39l*ll+-03i 

-0.14051-02 

0.93660-31 

0.24696-04 

0.19117-041 

0.19572-03 

0.89543-051 

300 

-0.51+121+-03 

0.20962-031 

-0.10824-02 

0.56307-031 

0.20621-04 

0.15615 -o 41 

0.12232-03 

0.52230-051 

1+00 

-0.1+3026-03 

0.1299^-031 

-0.86694-03 

0.38555-031 

0.16891-04 

0.13207-041 

0.85630-04 

0.33593-051 

500 

-0.351+81-03 

0.8891+2-Ol+i 

-0.71922-03 

0.28691-031 

0 .i 4 i 31-04 

0.11392-041 

0.63610-04 

0.22385-051 

600 

-0.30107-03 

0.65131+-01+1 

-0.61312-03 

0.22572-031 

0.12089-04 

0.99882-051 

0.48930-04 

0.14923-05 
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TABLE B.l (Continued) 


Reduced Frequency =0.50 Mechanical Admittance Magnitude M = 0.7 


y 

H^(ik) 

HgUk) 

T^(ik) 

Tgdk) 

200 

-0.10192-03 

0.l8915-05i 

-O.2U138-O3 

0.13861-031 

-0.32162-05 

0.50880-051 

0.85685-OU 

0.52259-051 

300 

-0.68039-OU 

-0.31067-05i 

-0.16778-03 

0.83710-OUi 

-0.21969-05 

0.36082-051 

O.53088-OU 

0.23681-051 

uoo 

-O.50906-0U 

-0.39600-051 

-0.12808-03 

0.59281-OUi 

-0.16755-05 

0.278U1-051 

O.37559-0U 

0.19525-051 

500 

-0.U0627-0U 

-0.39U38-051 

-O.IO3U5-O3 

O.U569U-OU1 

-0.13556-05 

0.22637-051 

O.28037-OU 

0.13289-051 

6oo 

-0.33788-0U 

-0.37lUU-05i 

-O.86727-0U 

0.37100-OUi 

-0.11388-05 

0.19063-051 

O.21568-OU 

0.88818-051 

Reduced Frequency = 0.90 

Mechanical Admittance Magnitude M = 0.7 

y 


Hg(ik) 

T^(lk) 

Tgdk) 

200 

-O.U5535-OU 

0.93735-06i 

-O.77OI3-OU 

0.20722-031 

0.66737-06 

-0.13689-051 

O.60316-OU 

0.32612-051 

300 

-O.29952-OU 

-0.29319-061 

-O.56665-OU 

0.13378-031 

O.U05U5-06 

-0.87057-061 

O.37769-OU 

O.I9IIU-O5I 

Uoo 

-O.2226O-0U 

-0.5578U-06i 

-O.UUU58-OU 

0.98321-OUi 

0.28758-06 

-0.63918-061 

0.26563-OU 

O.I278U-O5I 

500 

-0.17697-OU 

-0.60566-061 

-O.36U89-OU 

0.77597-OUi 

0.22175-06 

-0.50531-061 

0.19837-OU 

0.91510-061 

600 

-O.1U681-0U 

-0.59206-061 

-O.30912-OU 

0.6U0U3-0Ui 

O.I8OOU-06 

-0.U179U-06 i 

O.15079-OU 

0.58833-061 
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TABLE B.l (Continued) 


Reduced Frequency =0.98 Mechanical Admittance Magnitude M = 0.7 


y 


HgCik) 

T^(ik) 

Tgdk) 

200 

-O.26769-OU 

-0.85188-051 

-O.67352-OU 

0.2111U-031 

-0.19635-05 

0.11527-051 

-O.31220-OU 

-0.16099-021 

300 

-O.I7OIO-OU 

-0.61722-051 

-O.U8032-OU 

0.133U6-031 

-0.13505-05 

0.77U67-061 

O.II395-0U 

-0.1030U-021 

Uoo 

-O.I245I-OU 

-0.U7936-051 

-O.370U9-OU 

0.97UU8-0U1 

-0.10279-05 

0.58266-061 

O.I9553-OU 

-0.75713-031 

500 

-0.98151-05 

-0.39089-051 

-O.30095-OU 

0.76711-OUi 

-O.829U3-06 

0.U6677-061 

0.20667-OU 

-0.59827-031 

6oo 

-0.80992-05 

-0.32969-051 

-O.25321-OU 

0.632U2-0U1 

-0.69508-06 

0.38928-061 

O.19928-0U 

-0.U9UU6-031 

Reduced Frequency = 1.00 

Mechanical 

Admittance Msignltude M = 0.7 

y 

H^(ik) 

Hgdk) 

Tj^dk) 

Tgdk) 

200 

O.2327U-O2 

0.83762-021 

-O.U951U-02 

0.630U1-011 

-0.27985-OU 

O.I2U72-OIH 

0.65219-00 

0.56613-OOi 

300 

0.10921-02 

0.36632-021 

-O.I9U6U-O2 

0.U1722-011 

-O.III83--OU 

0.89U67-051 

O.U1736-OO 

0.43271-001 

Uoo 

0.62282-03 

0.20396-021 

-0.10303-02 

0.31036-011 

-O.583U6-05 

0.58528-051 

0.29174-00 

0.21856-001 

500 

O.UOI2O-O3 

0.21975-021 

-0.63796-03 

0.2U67U-011 

-0.35132-05 

0.U010U-05i 

0.21651-00 

0.14233-001 

600 

0.28001-03 

0.89789-031 

-O.U3U82-O3 

O.2OU65-011 

-0.23111-05 

0.28828-051 

0.16256-00 

0.99056-Oli 
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TABLE B.l (Continued) 


Reduced Frequency = 1.02 Mechanics^. Admittance Magnitude M = 0.7 


y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-O.2389U-OU 

-0.100T6-0Ui 

-O.60252-0U 

0.19517-031 

-O.1782U-05 

0.1531+5-051 

-0.28668-04 

-0.14991-021 

300 

-O.I5167-OU 

-0.70986-051 

-0.1*2758-01* 

0.12335-031 

-0.12266-05 

0.10246-051 

0.84268-05 

-0.95861-031 

4oo 

-0.11097-OU 

-0.51*506-051 

-0.32918-01* 

0.9007l*-0Ui 

-0.93381-06 

0.76857-061 

0.15748-04 

-0.70408-031 

500 

-O.87U57-O5 

-0.UUl7U-05i 

-0.26713-01* 

0.70911-01*1 

-0.75356-06 

0.6 i 479-06 i 

0.16906-04 

-0.55623-031 

600 

-0.72157-05 

-0.37115-051 

-O.22U6I-OI* 

0.581*62-01*1 

-0.63153-06 

0.51225-06 

0 .i 6409-04 

-0.45966-031 

Reduced Frequency = 1.10 

Mechanical 

Admittance Magnitude M = 0.7 

y 

H^(ik) 

Hgdk) 

T^(ik) 

Tgdk) 

200 

-0.31*650-03 

-0.59857-051 

-O.5278U-OI* 

0.29375-031 

0.23250-05 

-0.12770-061 

-0.50639-04 

-0.25669-051 

300 

-O.2221O-OU 

-0.1*281*7-051 

-0.381*63-01* 

0.18105-031 

0.15104-05 

-0.99077-071 

-0.31743-04 

-0.15174-051 

Uoo 

-0.1631*9-01* 

-0.3311*8-051 

-0.29938-01* 

0.3107l*-03i 

0.11189-05 

-0.78953-071 

-0.22631-05 

-0.99674-061 

500 

-0.12936-01* 

-0.26982-051 

-O.2UUU2-OI* 

0.10227-031 

0.88864-06 

-0.65264-071 

-0.17436-04 

-0.72613-061 

6oo 

-0.10702-0I* 

-0.22735-051 

-0.20630-01* 

0.83979-OUi 

0.73701-06 

-0.55514-071 

-0.12459-04 

-0.4 o 557-061 
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TABLE B.l (Continued) 


Reduced Frequency = 1.50 Mechanical Admittance Magnitude M = O.T 


y 

•H 

Hgdk) 

T^(ik) 

T^dk) 

200 

-0.17556-0U 

-0.1+8599-051 

0.16482-03 

0.45628-031 

0.49043-06 

-0.22569-051 

-0.35726-04 

-0.23157-051 

300 

-0.10I+672-0U 

-0.31937-05i 

0.80551-04 

0.26402-031 

0.29199-06 

-0.14516-051 

-0.22611-04 

-0.14556-051 

ItOO 

-0.71+787-05 

-0.23668-05i 

0.52141-04 

0.18556-031 

0.20928-06 

-0.10727-051 

-0.16378-04 

-0.98635-061 

500 

-0.58227-05 

-0.18781+-051 

0.38253-04 

0.14301-031 

0.16339-06 

-0.85132-061 

-0.12872-04 

-0.61467-061 

6oo 

-0.1+7688-05 

-0.15567-05i 

0 . 30107-04 

0.11633-031 

0.13410-06 

-0.70587-061 

-0.98889-05 

-0.42146-061 

Reduced Frequency = 1.60 

Mechanical 

A( 3 mittance Magnitude M = 0.7 

y 

H^(ik) 

H^dk) 

T^(ik) 

T2(lk) 

200 

-0.89327-05 

-0.30113-05i 

0.16532-04 

0.68915-041 

0.62887-05 

-0.41736-051 

-0.33829-04 

-0.21616-051 

300 

-0.62165-01+ 

-0.19287-051 

0.10278-03 

0.43882-041 

0.37168-05 

-0.25322-051 

-0.21364-05 

-0.12373-051 

Uoo 

-0.1+1+472-05 

-0.l4893-05i 

0.75617-05 

0.30887-041 

0.25417-05 

-0.19663-051 

-0.14472-04 

-0.74479-061 

500 

-0.37614-05 

-0.10836-051 

0.49211-05 

0.22314-051 

0.17746-05 

-0.14103-051 

-0.10436-04 

-0.58225-061 

6oo 

-0.31152-05 

-0.87316-061 

0.34364-05 

0.17225-041 

0.12564-05 

-0.92747-061 

-0.83456-05 
-0 .30262-061 
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TABLE B.i (Continued) 


Reduced Frequency =1.67 Mechanical Admittance Magnitude M = 0.7 


y 

H^(ik) 

H^dk) 

T^(ik) 

T^dk) 

200 

-0.84668-05 

-0.28446-05i 

0.16 i 47-04 

0.66839-041 

0.61174-05 

-0.40236-051 

-0.33226-04 

-0.19987-051 

300 

-0.60812-05 

-0.18981-051 

0.99746-05 

0.43l6l-04i 

0.36814-05 

-0.24132-051 

-0.21748-04 

-0.12126-051 

Uoo 

-0.42864-05 

-0.l4l75-05i 

0.74287-05 

0.30336-041 

0.24278-05 

-0.19113-051 

-0.16573-04 

-0.78451-061 

500 

-0.36214-05 

-0.10366-051 

0.48007-05 

0.22061 -o 41 

0.17542-05 

-0.31905-051 

-0.11768-04 
-0 . 54362-061 

600 

-0.29414-05 

-0.84i67-06i 

0.32874-05 

0.16932-041 

0.12217-05 

-0.91286-061 

-0.95359-05 

-0.34231-061 

Reduced 

Frequency = 2.00 

Mechanical 

Admittance Magnitude M = O.7 

y 

H^(ik) 

H2(ik) 

Tj^dk) 

Tgdk) 

200 

-0.86315-05 

-0.25004-05i 

0.16224-04 

0.57218-041 

0.12568-06 

0.69552-061 

-0.30172-04 

-0.18248-051 

300 

-0.51647-05 

-0.19237-051 

0.10117-04 

0.36238-041 

0.81715-07 

0.46158-061 

-0.18775-04 

-0.10823-051 

4oo 

-0.39713-05 

-0.13772-051 

0.81793-05 

0.27184-051 

0.56236-07 

0.36143-061 

-0 .i 4827-04 

-0.68338-061 

500 

-0.29117-05 

-0.10216-051 

0.63177-05 

0.19735-041 

0.43756-07 

0.28552-061 

-0.99754-05 

-0.42663-061 

600 

-0.22188-05 

-0.86504 -o 6 i 

0 . 52119-05 

0.14823-041 

0.35713-07 

0.23462-061 

-0.74611-05 

-0.29793-061 
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